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ABSTRACT 



An investigation of the interaction effects of a lower heated tube on pool boiling of 
pure R-124 from an upper horizontal tube was conducted at a saturation temperature of 
2.2 °C. The test tubes used were: (1) smooth tubes and (2) deformed surface (TURBO-B) 
enhanced tubes. The effects of tube spacing/configuration and lower tube heat flux on the 
heat transfer performance of the upper tube were investigated. 

For both tube arrays, the enhancing effect of bubbles from a lower tube was 
dramatic. This enhancement increased as lower tube heat flux increased. However, when 
upper tube heat fluxes were greater than 20 kW/m^ all enhancement disappeared. 

For a smooth tube array in natural convection, the effect of a lower heated tube on 
the heat transfer from an upper tube was small. In nucleate boiling, a P/D of 1.8 gave the 
best upper tube heat transfer performance and a vigorously nucleating lower tube 
eliminated upper tube hysteresis. With the lower tube unheated and an upper tube heat 
flux of greater than 3 kW/m^ the performance using R-124 was generally better than for 
R-1 14. With a nucleating lower tube (at 10 kW/m^), again the performance of R-124 was 
better, but only for upper tube heat fluxes of greater than 40 kW/m^. 

For a TURBO-B tube array, a 30 degree offset of the upper tube reduced the upper 
tube heat transfer performance (compared to the in-line configurations). This may 
indicate bubbles depart TURBO-B tubes differently than smooth tubes. 
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NOMENCLATURE 



This nomenclature reflects that which is found within 
programs SETUPJY (Appendix A) and DRPJY (Appendix B) . Note 



that 


the upper 


and lower test tube dimensions are the same 


since 


like test 


tubes are used. 
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I . INTRODUCTION 



The United States Navy intends to eliminate its dependence 
on chlorof luorocarbons (CFC's) as soon as feasible. 
Currently, the Navy is using CFC-114 (Refrigerant or R-114) in 
most of it's chilled water air conditioning (AC) plants. 
These shipboard AC plants provide the chilled cooling water to 
cool the vital weapons systems and electronic components 
besides providing habitability comfort for the crew. 

The Navy has identified a hydrochlorofluorocarbon, HCFC- 
124 or R-124, as a possible mid-term replacement for R-114 in 
naval shipboard AC systems. With an ozone depletion potential 
(ODP) of only 2% compared to 70% for R-114, R-124 is a much 
more tolerable fluid to use until an acceptable alternative 
refrigerant can be found that has zero ODP. 

Although both refrigerants operate near the same 
evaporating temperature (2.2 °C) , the characteristics of R-124 
under similar conditions as R-114 are quite different. For 
example, the refrigeration components (i.e., heat exchangers, 
piping, etc.) must be able to withstand: (1) the higher 
operating pressures of R-124 (11 psig at 2.2 °C compared to 0 
psig for R-114) and (2) the higher vapor pressure after the 
system is secured (in a main engineering space at 100 °F, R- 
124 would reach 66 psig compared to 31 psig for R-114) . Given 
the system components can operate safely at these higher 
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pressures of R-124, then the hardware modifications to the 
system should be limited to: 



1. Modifying the compressor to maintain rated capacity and 
use a miscible alkylbenzene oil for lubrication. 

2. Reset/replace the system safety devices to account for 
the higher pressures of R-124. 

3. Replace the system seals with neoprene seals, which 
resist deterioration when in contact with R-124, to 
maintain leak tight integrity. 



Besides hardware modifications, another main concern is 
how does the heat transfer characteristics of R-124 
compare to R-114. Table 1, from Bertsch [Ref. 1], provides 

TABLE 1. R-114/R-124 PHYSICAL PROPERTIES AT 2 . 2 °C 



Property 


R-114 


R-124 ' 


Psat (kPa) 


96 


177 


p(liq) (kg/mM 


1526 


1429 


p (vapor) (kg/m^) 


7.5 


11.3 


Cp(liq) (J/kg-K) 


932 


1060 


hfq (kJ/kg) 


134.7 


158.9 


k(liq) (W/m-K) 


0.0701 


0.0744 


p(liq) (g/m*s) 


0.4449 


0.3375 


0 (N/m) 


0.0136 


0.0129 



the physical properties of R-114 and R-124 at a saturation 
temperature of 2.2 °C. As seen from Table 1, R-124 has 

higher latent heat of vaporization, vapor density and 
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thermal conductivity than R-114 at 2.2 °C and a considerably 
lower liquid viscosity. Therefore, it is reasonable for one 
to assume that R-124 would have better heat transfer 
performance at 2.2 °C than R-114. Bertsch [Ref. 1] evaluated 
and compared the heat transfer performance of R-114 and R-124 
for a single horizontal smooth tube and various enhanced tubes 
using pure refrigerant and refrigerant/oil mixtures. One of 
his conclusions was: 

For smooth and finned tubes, significant improvements in 
heat transfer (>50%) were obtained in pure R-124 compared 
with pure R-114 at the same saturation temperature. This 
was attributed to the higher saturation pressure of R-124 
activating more nucleation sites and the higher latent 
heat of vaporization, vapor density and thermal 
conductivity of R-124 over R-114. 

Perry [Ref. 2] investigated the effect of oil, pressurization 
and subcooling on the onset of nucleate boiling of R-124 from 
a single horizontal tube. When Perry [Ref. 2] compared his 
results to Bertsch [Ref. 1] for pure R-124, he found a larger 
incipient wall superheat was being obtained, but almost 
identical data trends were observed at high heat fluxes and 
decreasing heat fluxes. He attributed the larger incipience 
to his degassing procedure of the refrigerent system after 
reassembly. 

Now, what if multiple tubes are arranged in a bundle 
configuration? How do the lower tubes influence the heat 
transfer performance of the upper tubes? Lake [Ref. 3] in 
search of an answer to these questions, used a configuration 
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of two horizontal tubes (smooth and enhanced tubes) in 
vertical alignment, submerged in a pool of R-114, separated by 
a pitch-to-diameter (P/D) ratio. After using various 

combinations of P/D ratios and lower tube heat fluxes. Lake 
[Ref. 3] found that a nucleating lower tube significantly 
enhanced the heat transfer performance of an upper tube, 
particularly at low heat fluxes. Also, a P/D ratio of 1.8 
gave the best heat transfer performance for smooth tubes at 
high heat fluxes (no P/D effect was found using the enhanced 
tubes) . 

To fully understand the nucleate pool boiling 

characteristics of R-124 and to directly compare the heat 
transfer performance of R-124 to R-114, multiple tube pool 
boiling data for R-124 is needed. Therefore, the following 
thesis objectives were established: 

1. Collect natural convection and pool boiling heat transfer 
data of pure R-124 using a single smooth and enhanced 
tube and check the repeatability with the existing R-124 
data of Bertsch [Ref. 1] and Perry [Ref. 2]. 

2. Collect natural convection and pool boiling heat transfer 
data of R-124 over a range of increasing and decreasing 
heat fluxes using two tubes while varying the spacing 
between the tubes and the heat flux on the lower tube. 

3. Compare the R-124 data to the R-114 data obtained by Lake 
[Ref. 3] . 
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II. LITERATURE SURVEY 



A. POOL BOILING 

Pool boiling occurs when a heated surface is submerged 
below the free surface of a liquid pool, at a temperature 
sufficiently above the local saturation temperature to cause 
bubbles to form from the surface. Figure 2.1, from Bar-Cohen 
[Ref. 4], is a pool boiling curve for a highly wetting liquid 
with several regimes displayed. Only the natural convection 
(or single phase convection) , boiling incipience and nucleate 
boiling regimes are discussed. 

1 . Natural Convection 

The regime (a-b) in Figure 2.1, is where heat transfer 
occurs because of the liquid flow caused by buoyancy forces 
due to density differences in the liquid pool. In this 
regime, the heat flux from the heated surface is proportional 
to some power of wall superheat (i.e., q" « or . 

This proportionality remains valid until the wall superheat 
reaches a value where vapor bubbles first form on the heated 
surface . 

2. Boiling Incipience 

The onset of nucleate boiling (ONB) or incipient 
boiling is a random phenomenon (You et al. [Ref. 5]) which 
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HEAT FLUX, q" (W/m^) 




WALL SUPERHEAT, 



Figure 2 . 1 Pool Boiling Curve 
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is influenced by factors such as surface condition of the 
heated surface, liquid properties and dissolved gases in the 
liquid pool. Bar-Cohen and Simon [Ref. 6] and Tong et al. 
[Ref. 7] investigated boiling incipience of highly-wetting 
liquids using refrigerants and dielectrics. They found that 
these liquids effectively flood the micro-cavities on a tube 
surface, causing fewer active embryos (micro-cavities 
containing trapped vapor) . Ultimately, these highly-wetting 
liquids require larger wall superheats to initialize bubble 
growth in the flooded micro-cavities. As a result of the 
larger wall superheats (up to 100 °C, Bar-Cohen and Simon 
[Ref. 6]), nucleation at (b) , in Figure 2.1, occurs in an 
explosive manner. This explosive boiling from the active 
micro-cavities causes the larger inactive cavities to 
activate, rapidly decreasing the required wall superheat to 
sustain boiling, noted as AT^^ in Figure 2.1. This drift from 
(b) to (c) only occurs during increasing heat flux and is 
referred to as thermal overshoot or superheat excursion. 

3 . Nucleate Boiling 

The regime (c-d) in Figure 2.1 is the fully developed 
nucleate boiling region. At (c) , small isolated bubbles are 
leaving the heated surface at low heat fluxes. As the heat 
flux is increased, these small bubbles become more numerous as 
more active sites are created. At high heat fluxes, these 
bubbles begin to coalesce, forming slugs of vapor, which grow 
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and depart from the surface. According to Thome [Ref. 8], in 
this regime, heat removal is accomplished by sensible heat in 
the form of superheated liquid or as latent heat in vapor 
bubbles. Also, the principal mechanisms of heat transfer are: 
(1) liquid-phase convection, enhanced by the stripping of the 
thermal boundary layer by the departing vapor bubbles, and (2) 
evaporation of the thin liquid micro-layer trapped between a 
growing vapor bubble and the heated surface. 



B. SINGLE TUBE HEAT TRANSFER EXPERIMENTS WITH R-124 AT NAVAL 

POSTGRADUATE SCHOOL 

Bertsch [Ref. 1] and Perry [Ref. 2] independently studied 
the heat transfer performance of single smooth and enhanced 
tubes submerged in a liquid pool of R-124 at 2.2 °C. Perry 
[Ref. 2] determined that in natural convection the 
experimentally based correlation, equation (1), developed by 
Churchill and Chu [Ref. 9] fairly well predicted the heat 
transfer performance of a smooth tube. 



NU = 



0.6 +0.387/ 



Ra 



\ (1 + (0.559/Pr)^/^®)^®/^ 



1/6 



12 



( 1 ) 



where Nu = Nusselt number = hDo/k 

Ra = Rayleigh number = (g'P-[Twall - Tsat } -Do^) / (u -a) 
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During nucleate boiling, most correlations available to 
predict the heat transfer performance of a smooth tube are 
expressed as heat flux as a function of tube wall superheat. 
Bertsch [Ref. 1] used a best fit correlation through his 
experimental nucleate boiling data to express the heat flux 
from a smooth single tube in pure R-124 as: 

q" = (17.0) AV (W/m^) (2) 

In the above equation, the wall superheat AT„ is expressed in 
degrees Celsius. 

For a single enhanced tube, Bertsch [Ref. 1] concluded 
that nucleate pool boiling enhancement (i.e., improvements in 
heat transfer compared to a single smooth tube) were as much 
as 6.2 times for re-entrant cavity surfaces (TURBO-B and HIGH 
FLUX) . Thome [Ref. 10] attributes these enhancements to the 
numerous re-entrant cavities that trap vapor, which create a 
large boiling site density, reducing the wall superheat 
required to sustain nucleation. 

C. SMALL TUBE ARRAY HEAT TRANSFER PERFORMANCE 

The effect of a lower tube on the heat transfer 
performance of an upper tube is very important when designing 
flooded evaporators. It is essential to understand the 
mechanisms that affect the heat transfer performance in small 
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tube arrays in natural convection and nucleate boiling for 
both smooth and enhanced surfaces. 

1 . Natural Convection 

In natural convection, the interaction effects of a 
lower heated tube on the heat transfer performance of an upper 
tube are the result of two offsetting mechanisms. With a 
lower tube heated, a rising heated buoyant plume impinges on 
and flows around an upper tube, which subjects the upper tube 
to local forced convection effects. This forced convection 
increases the local fluid velocity around the upper tube and 
raises the upper tube heat transfer coefficient. 
Additionally, the heated plume from the lower tube increases 
the local ambient temperature of the fluid impinging on the 
upper tube. This reduces the natural convection buoyancy 
driving force for this tube, decreasing the upper tube heat 
transfer coefficient. Therefore, the presence of a lower 
heated tube has opposing influences on the heat transfer 
performance of an upper tube. Which of these influences is 
dominant depends on the tube spacing, local heat flux and the 
thermophysical properties of the fluid. 

Sparrow and Niethammer [Ref. 11] studied the effect of 
vertical separation between two horizontal cylinders on the 
heat transfer performance of an upper cylinder in natural 
convection using air. It was discovered that heat transfer 
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enhancement^ increased as the distance between the tubes 
increased with a maximum enhancement of approximately 1.3 for 
7 < P/D < 9. Also, for a P/D < 3, the enhancement decreases 
below 1.0. Al-Alusi and Bushnell [Ref. 12] and Tokura et al. 
[Ref. 13] also working with a small array of horizontal smooth 
cylinders in air, found generally the same results. As the 
distance between the cylinders increased, the heat transfer 
enhancement increased, reaching a maximum, then returning to 
a lower value of enhancement. 

Sparrow and Boessneck [Ref. 14] investigated the 
effect of transverse misalignment of a lower cylinder on the 
heat transfer performance of an upper cylinder, for 2 ^ P/D ^ 
9, using air. For the perfectly aligned upper and lower 
cylinders, the heat transfer enhancement of the upper cylinder 
increased as the P/D ratio increased. Again, the largest 
enhancement was found with a P/D of 9 and the P/D of 2 
produced the lowest enhancment. As misalignment was 
increased, the enhancement began to merge to a common value of 
about 1.03 at sufficiently large offsets (about 2.5 cylinder 
diameters) . Sparrow and Boessneck concluded that this 
behavior is caused by the rising plume from the lower cylinder 
presenting a pre-heated quasi-forced-convection flow to the 



^Enhancement is defined for a small tube array as the 
upper cylinder heat transfer coefficient when both cylinders 
are heated divided by the upper cylinder heat transfer 
coefficient when the lower tube is unheated. 
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upper cylinder. At smaller P/D ratios, the increased local 
ambient temperatures of the fluid impinging on the upper 
cylinder dominates and decreases the upper cylinder heat 
transfer performance. At larger P/D ratios, the forced 
convection effects are more dominant and increase the heat 
transfer from the upper cylinder. 

Inagaki and Komori [Ref. 15] studied natural 
convection heat transfer from two vertically aligned 
horizontal cylinders in water. As discovered with air, a 
decrease in enhancement which falls below 1.0, is found for 
P/D < 1.0. A maximum enhancement of 1.25 occurs between 2.5 
^ P/D ^ 3.0 and at larger P/D ratios, the enhancement 
decreases to 1.0. Inagaki and Komori attribute this behavior 
to the same two opposing influences discussed previously. 

2 . Nucleate Boiling 

For small tube arrays in nucleate boiling, the most 
significant mechanism of enhancement within the array, is the 
well documented bundle effect^. Cornwell and associates 
[Refs. 16-19] studied fluid flow in tube bundles and found 
significiant enhancements in heat transfer in the upper 
portions of the bundles, concluding that forced convection and 
bubbly flow from the lower tubes, creating sliding bubbles 



^Bundle effect is the enhancement of an upper tube due to 
the activation (i.e., the nucleation of bubbles) of lower 
tubes within the bundle. 
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around the upper tube, cause significant enhancement at low 
heat flux conditions. 

Fujita et al. [Ref. 20] conducted tests with an array 
of three horizontal smooth tubes in R-113 at 0.1 MPa. For 
heat fluxes < 20 kW/m^, the bundle effect enhanced the upper 
tube heat transfer performance. Fujita et al. [Ref. 20] also 
used an array of two smooth tubes, where the lower tube heat 
flux was held constant and the upper tube heat flux was 
decreased uniformly. The upper tube enhancement increased as 
the lower tube heat flux was raised. Lake [Ref. 3] 
investigated the effect of P/D ratio and lower tube heat flux 
on the heat transfer performance of an upper tube, using an 
array of two smooth tubes in R-114 at 1 atmosphere. For these 
tests, the bottom tube heat flux was held constant and the 
upper tube heat flux was raised and lowered steadily. Lake 
found that for an upper tube heat flux of < 20 kW/m^, the 
enhancement of the upper tube was due to a bundle effect and 
the amount of enhancement increased with increased lower tube 
heat flux, which is consistent with Fujita et al . [Ref. 20]. 
Also, Lake found that for upper tube heat fluxes of > 20 
kW/m^, all bundle effects faded and that the P/D ratio of 1.8 
provided the best upper tube heat transfer performance. 

For a small array of enhanced tubes, several authors 
report findings similar to those found for the smooth tube 
arrays. Stephan and Mitrovic [Ref. 21], using a small array 
of three horizontal T-shaped finned tubes in R-114 at 0.2 MPa 
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and R-12 at 0.15 MPa, found bundle effects present at heat 
fluxes < 20 kW/m^, with these effects disappearing as the heat 
flux was raised. Hahne and associates [Refs. 22-24], using 
two different fin arrays in R-11 at 0.1 MPa, also found 
significant bundle effects for heat fluxes of < 20 kW/m^, 
which subside at higher heat fluxes. Hahne et al. [Ref. 24] 
also investigated the influence of P/D ratio on upper tube 
heat transfer. For heat fluxes < 20 kW/m^, the maximum P/D 
ratio tested (P/D of 3.0) produced the best heat transfer 
performance for both finned arrays. Above 20 kW/m^, neither 
P/D ratio nor tube fin pitch affected the upper tube heat 
transfer. Lake [Ref. 3], in conjunction with his testing of 
smooth tube performance above, found no effect of P/D ratio on 
HIGH FLUX tube (porous coated tube) performance at any heat 
flux and that there was no effect from a lower heated tube on 
the heat transfer performance of a nucleating upper HIGH FLUX 
tube. Fujita et al. [Ref. 20], like Lake, also tested porous 
coated tubes using the same testing conditions as for the 
smooth tubes and he concluded that the lower tube heat flux 
setting had no influence on the heat transfer coefficient of 
the upper tube, which agrees with the conclusions of Lake 
[Ref. 3]. 
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III. 



EXPERIMENTAL APPARATUS 



A. SYSTEM DESCRIPTION 

A schematic diagram of the experimental apparatus is shown 
in Figure 3.1. Details of the apparatus design, construction 
and calibration are discussed by Karasabun [Ref. 25], Reilley 
[Ref. 26] and Sugiyama [Ref. 27] . The apparatus is very 
similar to the one used by Lake [Ref. 3] to conduct pool 
boiling of pure R-114 with a two tube configuration. The only 
major change being the replacement of the Pyrex-glass Tees 
used as the condenser and evaporator. Bertsch [Ref. 1] 
replaced the Pyrex-glass Tees with fiberglass reinforced 
Pyrex-glass Tees enabling the system to withstand the 
increased operating pressure of R-124. 

As shown in Figure 3.1, the present system contains the 
following components: (1) a reinforced Pyrex-glass Tee, used 
as the evaporator to boil the refrigerant; (2) a reinforced 
Pyrex-glass Tee, used as the condenser to condense the 
refrigerant vapor; (3) a liquid refrigerant storage 
reservoir; (4) an oil storage and mixing subsystem (not 
used); (5) a vacuum pump; (6) an ethylene-glycol/water 
mixture cooling subsystem; (7) a condenser pump; 
(8) aluminum framework with plexiglass siding which houses 
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Figure 3.1 Scheiuatic of Experimental Apparatus 
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components (1) through (4); and (9) a data acquisition and 
instrumentation system (not shown in Figure 3.1). 

The system operates in a closed loop with vapor being 
generated by the heated tubes submerged in the R-124 liquid 
pool in the evaporator. The vapor travels to the condenser, 
via a two inch diameter aluminum L-shaped pipe, where it's 
condensed on a copper coil cooled by the chilled ethylene- 
glycol/water mixture. The condensed liquid refrigerant is 
gravity fed back to the evaporator via the return line. 

B. BOILING TEST SECTION 

1 . Evaporator 

The evaporator is a Corning Pyrex-glass Tee section 
(three inch nominal inner diameter) coated with a continuous 
winding of fiberglass impregnated filaments bonded with a 
modified epoxy resin. It was installed horizontally with the 
side-arm of the Tee vertical and hydrostatically safety tested 
to 75 psig by Bertsch [Ref. 1]. Aluminum endplates, 210 mm 
diameter and 12.7 mm thick, house the thermocouple wells, an 
oil entry fitting (not used) and the Teflon inserts, which 
hold the test tubes in position. The endplates are mated to 
the evaporator using cast-iron flanges. At the very bottom of 
the evaporator sits the refrigerant liquid return line from 
the condenser. A copper deflector plate was installed just 
above the return line to prevent any vapor bubbles that might 
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rise off the return line from impinging on the test tubes. The 
liquid refrigerant level in the evaporator is approximately 10 
mm below the top of the evaporator and approximately 10 mm 
above the top of the upper tube. Figure 3.2 is a sketch of 
the evaporator and a cast-iron flange. 

2 . Evaporator Tubes 

Figure 3.3 is a schematic of a boiling tube. Two 
different tube types were used in this study: a smooth tube 
and a TURBO-B enhanced surface tube. Table 2 provides the 
dimensions of each tube type while Figure 3.4 shows an 



TABLE 2. DIMENSIONS OF TEST TUBES 



Tube 


D1 


D2 


Di 


Do 


L 


Lu 


kc 




mm 


mm 


mm 


mm 


mm 


mm 


W/m*K 


smooth 


12.4 


15.9 


12.7 


15.9 


203.2 


76.2 


344 


TURBO-B 


11.6 


13.8 


11.8 


15.9 


203.2 


‘ 76.2 


398 



expanded view of the enhanced surface of the TURBO-B tube. 
Each tube is heated by a 1 kW, 240 volt stainless-steel 
cartridge heater. The heater measures 6.35 mm outside 
diameter, 203.2 mm overall length with a heated length of 190 
mm and is centered inside the test tube. The tubes are 
instrumented with eight thermocouples at various locations to 
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Figure 3.2 Sketch of Pyrex Glass Evaporator 
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Figure 3.3 Schematic of a Test Tube 






Figure 3 . 4 Surface Profile of Enhanced (TURBO-B) Tube 
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monitor the tube wall temperatures. Bertsch [Ref. 1] provides 
a detailed assembly procedure for the manufacturing and 
instrumentation of a test tube. 

The test tubes are secured in the evaporator by Teflon 
inserts and sealed with o-rings (Figure 3.5). Strong-baclcs, 
shown in Figure 3.6, are used to anchor the Teflon inserts in 
the endplates. 

3 . Tube Configurations 

Figure 3.7 shows the six tube configurations that were 
used. Configuration CO is the single tube configuration used 
by Bertsch [Ref. 1] and Perry [Ref. 2] . In configurations Cl 
through C4, the top tube is fixed relative to the evaporator 
while the lower tube position is varied. Configuration C5 
consists of rotating the Teflon insert used for C4 through an 
angle of 30°. Note that when rotation is completed, the upper 
tube position has moved off the evaporator centerline and the 
lower tube position remains unchanged. 

C . CONDENSER SECTION 

The condenser is also a fiberglass-epoxy coated Corning 
Pyrex-glass Tee section identical in size, shape and substance 
to the evaporator shown in Figure 3.2. It is mounted 
vertically (see Figure 3.1) with aluminum endplates bolted to 
cast-iron flanges, similar to the evaporator arrangement. The 
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Figure 3 . 5 Sketch of an Endplate and Teflon Insert 
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Figure 3.6 Sketch of a Strongback 
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TUBE; O - NO TUBE (PLUGGED) 



top endplate is fitted with connections for refrigerant 
charging/vacuum operations and the ethylene-glycol/water 
coolant inlet tube. The bottom endplate is fitted with 
connections for the gravity feed return line to the evaporator 
and the ethylene-glycol/water coolant outlet tube. The 
refrigerant vapor rises from the evaporator, passing through 
the L-shaped aluminum pipe and enters the condenser where it 
is condensed on a 3/8 inch outside diameter copper tube which 
is helically coiled inside the condenser. 

D. REFRIGERANT STORAGE RESERVOIR 

The refrigerant storage reservoir is an aluminum cylinder, 
230 mm in diameter and 254 mm in height. It is located within 
the apparatus enclosure above the evaporator and below the 
condenser (see Figure 3.1). The reservoir enables the 
operator to store the liquid R-124 instead of releasing it 
into the atmosphere, prior to changing out the test tubes. 
The installed sight glass allows monitoring of the liquid 
level during replenishing from an outside source (i.e., 145 
lb. storage cylinder) or refilling the evaporator after a test 
tube changeout. The reservoir is isolated from or aligned to 
the system by means of valves V6, V7, V12 and V13. 
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E. OIL SUBSECTION 



This study was devoted to pure R-124. In order to ensure 
that no oil contaminated the system, the oil line from the oil 
cylinder to the evaporator was removed (see Figure 3.1) and 
the access on the evaporator endplate was capped. Perry 
[Ref. 2] describes the oil subsection components and the oil 
addition procedure for a refrigerant/oil mixture study. 

F . COOLING SUBSYSTEM 

1 . Ethylene-Glycol/Water Coolant Sump 

At the base of the apparatus sits a 0.15 cubic meter 
plexiglass tank containing approximately 30 gallons of 
ethylene-glycol/water mixture producing a freezing point of 
-25 °C. The sump rests on a wood platform and is wrapped with 
22 mm thick foam rubber sheet for insulation. Supply and 
discharge connections for the condenser and sump pump as well 
as a thermocouple to monitor the sump temperature are 
installed. 

2 . Refrigeration Plants 

A 1/2 ton R-502 and 1/4 ton R-12 refrigeration plants 
are used simultaneously to maintain the ethylene-glycol/water 
coolant temperature at about -15° C. Both plants are located 
close to the sump and contain the following components: (1) 
hermetically sealed compressor assembly; (2) air-cooled 
condenser; (3) receiver; (4) filter-dryer; (5) pressure 
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regulator; (6) temperature control switch; and (7) 
thermostatic expansion valve. The R-502 plant (Figure 3.8) 
uses a counter flow heat exchanger as the system evaporator 
while the R-12 plant uses a 9.5 mm diameter coiled copper 
tubing submerged in the sump tank as the system evaporator 
( Figure 3.1). 

3 . P\imps and Control Valve (VC) 

A positive displacement sump pump, rated at 8 gpm, is 
connected to the discharge connection of the coolant sump with 
a one inch diameter PVC pipe. It circulates coolant from the 
sump through the R-502 evaporator and back to the sump. A 
similar pump, called the condenser pump, circulates coolant 
from the sump to the condenser, via the control valve VC, and 
back to the sump (see Figure 3.1) . Valve VC is used by the 
operator to control the refrigerant pool saturation 
temperature. A manually operated bypass valve (V9, Figure 
3.1) is installed upstream of valve VC to prevent overloading 
the condenser pump when valve VC is closed. 



G. APPARATUS ENCLOSURE 

A welded aluminum frame measuring 214 cm high, 61 cm wide 
and 51 cm deep is used to contain most of the system 
components. The upper half of the frame contains the 
evaporator, condenser, oil subsystem, and refrigerant storage 
reservoir (see Figure 3.1). The lower half of the frame 
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Figure 3.8 Schematic of R-502 Refrigeration System 



contains the coolant sump tank. For operator safety and to 
reduce heat leakage into the evaporator and condenser, all 
sides of the upper half of the frame are covered with 13 mm 
thick plexiglass, hinged on both sides to allow access to the 
system components. Holes in the front plexiglass cover allow 
operation of valves VI through V8 while maintaining the 
integrity of the enclosed upper frame. 

H . INSTRUMENTATION 

1 . Power Measurement System 

Figure 3.9 is a schematic of the power measurement 
system. Each boiling tube heater is powered by a 240 volt AC 
source controlled by a variac. The variac output ranges from 
0-220 volt and 0-5 amp, adjustable by the operator to obtain 
the desired heat flux at the tube surface. The current and 
voltage sensors convert the AC input signal into DC voltages 
which are used as input by the Data Acquisition Unit. 

2 . Temperature Measurement 

Temperatures are measured using 30 gauge (0.25 mm 
diameter) copper-constantan (type T) thermocouples. For each 
test tube, eight thermocouples are spaced over the heater 
sleeve, as shown in Figure 3.10, to measure the wall 
temperatures. Three thermocouples are inserted into housings 
that penetrate the evaporator endplate and are submerged in 
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SCCJlON A-A (DIMENSIONS IN MM I 



the refrigerant liquid pool at different locations (see Figure 
3.1). Each housing is constructed of a low conductivity 
stainless steel body and a high conductivity copper tip 
(Figure 3.11). This minimizes the errors in pool temperature 
measurements caused by heat conduction through the housing. 
A single thermocouple is located in the coolant sump tan)c to 
measure the sump temperature. 

Thermocouples were read by the a Hewlett-Packard (HP) 
3497A Data Acquisition Unit which is controlled by an HP 326 
Computer System. 

I. DATA ACQUISITION AND REDUCTION PROGRAM 

Data acquisition and reduction program DRP72, used by Lake 
[Ref. 3], was modified to incorporate the thermophysical 
properties of saturated R-124. Bertsch [Ref. 1] discusses the 
generation of the equations used in the program calculations. 
The modified DRP72 program was renamed DRPJY. Appendix B 
contains the program DRPJY. The program DRPJY directs the HP 
362 Computer System and the HP 3497A Data Acquisition Unit to 
read and process all thermocouple, current and R.M.S. voltage 
readings. Using this information, DRPJY performs the 
following calculations (from Bertsch [Ref. 1]): 

1 . Converts voltage readings into tube temperatures and 
power . 

2. Computes heat transfer rate from cartridge heaters. 
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Figure 3.11 Sketch of Thermocouple Housing 



3. Computes average test tube wall temperatures and wall 
superheat (Twall - Tsat) . 

4. Computes physical properties of R-124 using property 
correlations from RefProp [Ref. 28] at the film 
temperature (Twall + Tsat)/2. 

5. Computes heat transfer coefficient of R-124 from the 
unheated ends and heated length of the test tubes. 

6. Computes heat losses from the unheated ends of the test 
tubes . 

7. Computes heat flux from heated length of test tubes. 



The processed output data was printed by an HP Inkjet 
printer and stored on the computer hard drive. As a 
precaution, each output data file was copied to a 3.5 inch 
floppy disk in the event the computer hard drive becomes 
unreadable. Table 3 lists the channel assignments that are 
sensor inputs to the Data Acquisition Unit. 
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TABLE 3. HP-3497A CHANNEL ASSIGNMENTS 



Channel 


Channel Assignments 


00-07 


Upper test tube wall temperature (T1-T8) 


08-15 


Lower test tube wall temperature (T9-T16) 


16 


Liquid pool thermocouple (Tldl) 


17 


Liquid pool thermocouple (Tld2) 


16 


Liquid pool thermocouple (Tv) 


19 


Sump coolant temperature 


2n 


R.M.S. voltage of tube heaters 


21 


Upper tube current sensor 


22 


Lower tube current sensor 
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IV. 



EXPERIMENTAL PROCEDURE 



A. EVAPORATOR DISASSEMBLY PROCEDURE 

The following procedure is used to transfer the liquid 
refrigerant from the evaporator to the reservoir to facilitate 
test tube changeouts: 



1. Operate both R-502 and R-12 refrigeration systems to 
cool the sump tank temperature to approximately -15 °C. 

2. Start the condenser pump and using control valve VC, 
adjust the coolant supply pressure to 1 psig. 

3. Close valve V5 to isolate the evaporator from the 
condensing refrigerant vapor (see Figure 3.1). 

4. Align the refrigerant storage reservoir to the system by 
opening valves V6, V7, V12 and V13. 

5. Position and energize the high volume air dryer (set for 
high heat) so the heated air stream is blowing on the 
center of the evaporator. Adjust valve VC to ensure the 
refrigerant vapor is condensing. 

6. After 1-2 hours, the evaporator should be dry with the 
liquid refrigerant in the reservoir and a small amount 
in the return line. 

7. Close valves V6, V7, V12 and V13, isolating the 
reservoir from the system. 

8. Secure and remove the air dryer. 

9. Stop the condenser pump. 

10. Open valves V8 and VIO to reduce to system pressure to 
atmospheric. Close valve V8 and VIO when the system 
pressure reads 0 psig on the system pressure gage. 

11. Disconnect the test tube heater and thermocouple leads. 
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12. Remove strongbacks. 

13. Remove Teflon inserts and test tubes. 

B. EVAPORATOR ASSEMBLY PROCEDURE 

1 . Evaporator Assembly 

Before installing smooth test tubes, rinse with 
acetone then wipe clean with lint-free wipes or if installing 
TURBO-B test tubes, brush the enhanced surfaces with a soft 
bristle brush (i.e., toothbrush) using ethyl alcohol then 
repeat the brushing using acetone and let the tubes air dry. 
Check the Teflon insert o-rings for damage and replace as 
necessary. Before installing the test tubes, wipe the 
evaporator interior clean with lint-free wipes. Install the 
test tubes in the evaporator using the Teflon inserts to 
secure the tubes in place. Snug the strongbacks up to the 
Teflon inserts and connect the thermocouple leads. Prior to 
connecting the heater leads, check the resistance between the 
heater leads and the outer tube wall using a multimeter. One 
megohm resistance or greater is satisfactory. 

2 . Leak Test 

After assembling the evaporator and before charging 
the system with refrigerant, perform a vacuum test to verify 
the integrity of the system. Evacuate the system by starting 
the vacuum pump and opening valves VI 1 and V8 (see Figure 
3.1). After evacuating the system to approximately 29 
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in. Hg, close valves Vll and V8, then stop the vacuum pump. 
If no drop in vacuum is detected (i.e., no leak) within 30 
minutes, charge the system with R-124. 

If a leak is detected, open valves V8 and Vll and 
pressurize the system with air to 20 psig. Apply a soap/water 
solution to the system sealing surfaces to identify the leak 
location (bubbles will be seen from the leak source) . After 
identifying and repairing the leak, retest the system for 
leaks . 

3. Charging System with Refrigerant 

Once the system integrity is assured, fill the 
evaporator with R-124, using one of the following procedures, 
to a level just helov the top of the evaporator ( NOTE : some 
refrigerant will be lost during the degas procedure, 
therefore, it is important to fill the evaporator above the 
level needed to conduct data runs) : 
a. From Storage Reservoir 

Open valves V12 and V7 to equalize the pressure 
between the storage reservoir and the system. Open valves 
V13, V6 and V5 to gravity fill the evaporator with liquid 
refrigerant. When the evaporator is filled to the proper 
level, close valves V6, V13, V7 and V12 to isolate the storage 
reservoir from the system. 
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b. From Refrigerant Cylinder 

Connect the cylinder to valve VI 0 using the 
charging hose (Figure 3.1). Open valves V5, V8 and VIO and 
then open the cylinder liquid valve to release the refrigerant 
from the cylinder into the system. Close valves V8, VIO and 
the cylinder liquid valve when the evaporator is filled to the 
proper level. Disconnect the charging hose from valve VIO and 
stow the cylinder. 

4 . Degassing Procedure 

This degas procedure is similar to Perry's [Ref. 3] 
degassing procedure B. Degas the system by boiling the 
refrigerant, using a upper test tube heat flux of 85 KW/m^ and 
a lower test tube heat flux of 20 KW/m^, for 30 minutes. 
Ensure the condenser pump is running and control valve VC is 
open to prevent overpressurizing the system. During the 
boiling process, air and non-condensible gases rise in the 
system and collect in the condenser. After the 30 minute 
boiling time has elapsed, secure the heat flux to the test 
tubes and stop the condenser pump. Let the system set for 30 
minutes allowing the condensed refrigerant to drain to the 
evaporator. After the 30 minute set time has elapsed, start 
the vacuum pump and open valve VI 1, then open valve V8 for 20 
seconds to remove the accumulated air and non-condensible 
gases. After closing valve V8, stop the vacuum pump and close 
valve VI 1. 



40 



5. Data Acquisition Channel Check 

After all evaporator thermocouples (test tubes and 
pool) are connected, load and run program SETUPJY to verify 
the output of each thermocouple is correct. Any output that 
is inconsistent or suspect is either fixed or noted and later 
excluded from the calculations in program DRPJY. 



C. STANDARD OPERATING PROCEDURE 

The following procedure is used to conduct each data 
collection run: 



1. Operate the R-502 and R-12 refrigeration systems to cool 
the sump to approximately -15 °C. 

2. Energize the HP-362 Computer System, Data Acquisition 
Unit and the variac control panel. 

3. Load and run program SETUPJY. Verify the thermocouple 
outputs are consistent. 

4. Verify control valve VC is closed, then start the 
condenser pump. 

5. Slowly open valve VC until the coolant pressure gage 
reads 1/8 psig. This pressure setting will ensure that 
the refrigerant pool will be slowly cooled. 

6. Monitor the thermocouple outputs with SETUPJY until 
confident the thermocouples are operating properly. 
Note which thermocouples are inconsistent and must be 
excluded from program DRPJY. 

7. When the refrigerant pool temperature is approaching 5 
°C, load and run program DRPJY. Select 'TAKE DATA' 
option . 

8. At the prompt, enter the filename that the data is to be 
stored under (page 38 explains the file naming scheme 
used) . 
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9. At the prompt, enter the number of defective upper tube 
thermocouples that were identified by SETUP JY. 

10. At the prompt, enter the defective upper tube 
thermocouple locations that were identified by SETUPJY. 

11. At the prompt, enter the number of defective lower tube 
thermocouples that were identified by SETUPJY. 

12. At the prompt, enter the defective lower tube 
thermocouple locations that were identified by SETUPJY. 

13. Select the type of test tubes installed. 

14. From the main menu, select 'SET Tsat' . Monitor the pool 
thermocouples and, if required, adjusted valve VC to 
maintain the desired pool saturation temperature of 2.22 
± 0.1 °C. 

15. From the main menu, select 'SET UPPER TUBE HEAT FLUX', 
then adjust the main variac to the desired heat flux. 
Initial heat flux setting is 500-600 KW/m^. 

16. Press the keyboard button fl to return to the main menu 
and select 'SET LOWER TUBE HEAT FLUX'. Adjust the aux 
variac to the desired heat flux. 

17. From the main menu, select 'SET Tsat'. Adjust valve VC 
to maintain the desired pool saturation temperature 
while monitoring the pool thermocouples. 

18. Once the pool saturation temperature is stable, return 
to the main menu and select 'TAKE DATA'. 

19. For each heat flux setting steps 15-18 are repeated. 



D . DATA FILES 

Data were filed using the following convention: 

Example (D0904ASC1) : 

1. First letter D indicates the file is a data file. 

2. The following four characters indicate the date; in this 
example September 04, 1993. 
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3. The next letter indicates the order in which the data run 
was completed on that day; no letter - 1st run, A - 2nd 
run, B - 3rd run, etc. 

4. The next one or two letters indicate the tubes used; 

S - smooth, TB - TURBO-B. 

5. The last letter and number combination represent the tube 
configuration used. 



Table 4 is a summarized list of the data runs completed. 
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TABLE 4. LIST OF DATA RUNS 



Data File 


Tube Type 


Purpose 


D0728ASC0 


Smooth 


Repeatability; single tube 


D0804ASC0 


Smooth 


Repeatability; single tube 


D0816SC0 


Smooth 


Repeatability; single tube 


D0520SC1 


Smooth 


Repeatability; P/D = 2.0; no 
power on lower tube 


D0714SC1 


Smooth 


Repeatability; P/D = 2.0; no 
power on lower tube 


D0831SC1 


Smooth 


Data; repeatability; P/D = 2.0; 
no power on lower tube 


D0901SC1 


Smooth 


Data; P/D =2.0; 1 kW/m^ on 
lower tube 


D0902SC1 


Smooth 


Data; P/D = 2.0; 5 kW/m^ on 
lower tube 


D0902BSC1 


Smooth 


Data; P/D = 2.0; 20 kW/m^ on 
lower tube 


D0903SC1 


Smooth 


Data; P/D = 2.0; q" on both 
tubes inc/dec together 


D0904SC1 


Smooth 


Data; P/D = 2.0; 5 kW/m^ on 
upper tube 


D0904ASC1 


Smooth 


Data; P/D = 2.0; 10 kW/m^ on 
upper tube 


D0905SC1 


Smooth 


Data; P/D = 2.0; 20 kW/m^ on 
upper tube 


D0906SC2 


Smooth 


Data; P/D = 1.8; no power on 
lower tube 


D0906ASC2 


Smooth 


Data; P/D = 1.8; 1 kW/m^ on 
lower tube 


D0907SC2 


Smooth 


Data; P/D = 1.8; 5 kW/m^ on 
lower tube 


D0907ASC2 


Smooth 


Data; P/D = 1.8; 10 kW/m^ on 
lower tube 
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TABLE 4. LIST OF DATA RUNS (continued) 



1 Data File 


Tube Type 


Purpose 


1 D0908SC2 


Smooth 


Data; P/D = 1.8; 20 kW/m^ on 
lower tube 1 


1 D0908ASC2 


Smooth 


Data; P/D = 1.8; q" on both 
tubes inc/dec together 


D0908ASC2 


Smooth 


Data; P/D = 1.8; 5 kW/m^ on 
upper tube 


D0909SC2 


Smooth 


Data; P/D = 1.8; 10 kW/m^ on 
upper tube 


D0909ASC2 


Smooth 


Data; P/D = 1.8; 20 kW/m^ on 
lower tube 


D0910SC3 


Smooth 


Data; P/D = 1.5; no power on 
lower tube 


D0910ASC3 


Smooth 


Data; P/D = 1.5; 1 kW/m^ on 
lower tube 


D0917SC3 


Smooth 


Data; P/D = 1.5; 5 kW/m^ on 
lower tube 


D0919SC3 


Smooth 


Data; P/D = 1.5; 10 kW/m^ on 
lower tube 


D0919ASC3 


Smooth 


Data; P/D = 1.8; 20 kW/m^ on 
lower tube 


D0920SC3 


Smooth 


Data; P/D = 1.8; q" on both 
tubes inc/dec together 


D0920ASC3 


Smooth 


Data; P/D = 1.5; 5 kW/m^ on 
upper tube 


D0920BSC3 


Smooth 


Data; P/D = 1.5; 10 kW/m^ on 
upper tube 


D0920CSC3 


Smooth 


Data; P/D = 1.5; 20 kW/m^ on 
upper tube 


D0921SC4 


Smooth 


Data; P/D = 1.2; no power on 
lower tube 


D0921ASC4 


Smooth 


Data; P/D = 1.2; 1 kW/m^ on 
lower tube 


D0921BSC4 


Smooth 


Data; P/D = 1.2; 5 kW/m^ on 
lower tube 
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TABLE 4. LIST OF DATA RUNS (continued) 



Data File 


Tube Type 


Purpose 


D0921CSC4 


Smooth 


Data; P/D = 1.2; 10 kW/m^ on 
lower tube 


D0922SC4 


Smooth 


Data; P/D = 1.2; 20 kW/m^ on 
lower tube 


D0922ASC4 


Smooth 


Data; P/D = 1.2; q" on both 
tubes inc/dec together 


D0922BSC4 


Smooth 


Data; P/D = 1.2; 5 kW/m^ on 
upper tube 


D0922CSC4 


Smooth 


Data; P/D = 1.2; 10 kW/m^ on 
upper tube 


D0922DSC4 


Smooth 


Data; P/D = 1.2; 20 kW/m^ on 
upper tube 


D0930SC5 


Smooth 


Data; P/D = 1.2 rotated 30°; no 
power on lower tube 


D0930ASC5 


Smooth 


Data; P/D = 1.2 rotated 30°; 1 
kW/m^ on lower tube 


D0930BSC5 


Smooth 


Data; P/D = 1.2 rotated 30°; 5 
kW/m^ on lower tube 


D0930CSC5 


Smooth 


Data; P/D = 1.2 rotated 30°; 10 
kW/m^ on lower tube 


D0929CSC5 


Smooth 


Data; P/D =1.2 rotated 30°; 20 
kW/m^ on lower tube 


D0929DSC5 


Smooth 


Data; P/D = 1.2 rotated 30°; q" 
on both tubes inc/dec together 


D0929BSC5 


Smooth 


Data; P/D = 1.2 rotated 30°; 5 
kW/m^ on upper tube 


D0929ASC5 


Smooth 


Data; P/D =1.2 rotated 30°; 10 
kW/m^ on upper tube 


D0929SC5 


Smooth 


Data; P/D = 1.2 rotated 30°; 20 
kW/m^ on upper tube 


D1012TBC0 


TURBO-B 


Repeatability; single tube 


D1012ATBC0 


TURBO-B 


Repeatability; single tube 


D1012BTBC0 


TURBO-B 


Repeatability; single tube 
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TABLE 4. LIST OF DATA RUNS (continued) 



Data File 


Tube Type 


Purpose 


D1019TBC1 


TURBO-B 


Repeatability; P/D = 2.0; no 
power on lower tube 


1 D1018ATBC1 


TURBO-B 


Repeatability; P/D = 2.0; no 
power on lower tube 


D1018BTBC1 


TURBO-B 


Data; repeatability; P/D = 2.0; 
no power on lower tube 


D1018CTBC1 


TURBO-B 


Data; P/D =2.0; 1 kW/m^ on 
lower tube 


D1018DTBC1 


TURBO-B 


Data; P/D =2.0; 5 kW/m^ on 
lower tube 


D1019TBC1 


TURBO-B 


Data; P/D =2.0; 10 kW/m^ on 
lower tube 


D1019ATBC1 


TURBO-B 


Data; P/D = 2.0; 20 kW/m^ on 
lower tube 


D1019ITBC1 


TURBO-B 


Data; P/D = 2.0; q" on both 
tubes decreased together 


D1019CTBC1 


TURBO-B 


Data; P/D =2.0; 5 kW/ra^ on 
upper tube 


D1019DTBC1 


TURBO-B 


Data; P/D = 2.0; 10 kW/m^ on 
upper tube 


D1019ETBC1 


TURBO-B 


Data; P/D = 2.0; 20 kW/m^ on 
upper tube 


D1020TBC2 


TURBO-B 


Data; P/D = 1.8; no power on 
lower tube 


D1020ATBC2 


TURBO-B 


Data; P/D =1.8; 1 kW/m^ on 
lower tube 


D1020BTBC2 


TURBO-B 


Data; P/D = 1.8; 5 kW/m^ on 
lower tube 


D1020CTBC2 


TURBO-B 


Data; P/D = 1.8; 10 kW/m^ on 
lower tube 


D1020DTBC2 


TURBO-B 


Data; P/D = 1.8; 20 kW/m^ on 
lower tube 
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TABLE 4. LIST OF DATA RUNS (continued) 



Data File 


Tube Type 


Purpose 


D1020ETBC2 


TURBO-B 


Data; P/D = 1.8; q" on both 
tubes decreased together 


D1021TBC2 


TURBO-B 


Data; P/D = 1.8; 5 kW/m^ on 
upper tube 


D1021ATBC2 


TURBO-B 


Data; P/D = 1.8; 10 kW/m^ on 
upper tube 


D1021ATBC2 


TURBO-B 


Data; P/D = 1.8; 20 kW/m^ on 
upper tube 


D1021CTBC3 


TURBO-B 


Data; P/D = 1.5; no power on 
lower tube 


D10211TBC3 


TURBO-B 


Data; P/D = 1.5; 1 kW/m^ on 
lower tube 


D1021ETBC3 


TURBO-B 


Data; P/D = 1.5; 5 kW/m^ on 
lower tube 


D1021FTBC3 


TURBO-B 


Data; P/D = 1.5; 10 kW/m^ on 
lower tube 


D1022TBC3 


TURBO-B 


Data; P/D = 1.5; 20 kW/m^ on 
lower tube 


D1022ATBC3 


TURBO-B 


Data; P/D = 1.5; q" on both 
tubes decreased together 


D1022BTBC3 


TURBO-B 


Data; P/D = 1.5; 5 kW/m^ on 
upper tube 


D1022CTBC3 


TURBO-B 


Data; P/D = 1.5; 10 kW/m^ on 
upper tube 


D1022DTBC3 


TURBO-B 


Data; P/D = 1.5; 20 kW/ra^ on 
upper tube 


D1022ETBC4 


TURBO-B 


Data; P/D = 1.2; no power on 
lower tube 


D1022FTBC4 


TURBO-B 


Data; P/D = 1.2; 1 kW/m^ on 
lower tube 


D1025TBC4 


TURBO-B 


Data; P/D = 1.2; 5 kW/m^ on 
lower tube 


D1025ATBC4 


TURBO-B 


Data; P/D = 1.2; 10 kW/m^ on 
lower tube 
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TABLE 4. LIST OF DATA RUNS (continued) 



i 

Data File 


Tube Type 


Purpose 1 


D1025BTBC4 


TURBO-B 


Data; P/D = 1.2; 20 kW/m^ on 
lower tube 


D1025CTBC4 


TURBO-B 


Data; P/D = 1.2; q" on both 
tubes decreased together 


D1025DTBC4 


TURBO-B 


Data; P/D = 1.2; 5 kW/m^ on 
upper tube 


D1025ETBC4 


TURBO-B 


Data; P/D = 1.2; 10 kW/m^ on 
upper tube 


D1025FTBC4 


TURBO-B 


Data; P/D = 1.2; 20 kW/m^ on 
upper tube 


D1026TBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; no 
power on lower tube 


D1026ATBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 1 
kW/m^ on lower tube 


D1026BTBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 5 
kW/m^ on lower tube 


D1026CTBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 10 
kW/m^ on lower tube 


D1026DTBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 20 
kW/m^ on lower tube 


D1026ETBC5 


TURBO-B 


Data; P/D =1.2 rotated 30°; q" 
on both tubes decreased 
together 


D1027TBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 5 
kW/m^ on upper tube 


D1027ATBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 10 
kW/m^ on upper tube 


D1027BTBC5 


TURBO-B 


Data; P/D = 1.2 rotated 30°; 20 
kW/m^ on upper tube 
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V. 



DISCUSSION OF RESULTS 



A. REPRODUCIBILITY/REPEATABILITY 

To verify that the equipment and operating procedure being 
used were consistent with those used by Bertsch [Ref. 1] and 
Perry [Ref. 2 ], six reproducibility/repeatability runs were 
performed. Program DRPGP (from Perry's single tube work 
[Ref. 2]) was used for the runs conducted with a single test 
tube and DRPJY was used for all runs conducted with two test 
tubes. Figures 5.1 and 5.2 are plots of increasing and 
decreasing heat flux using a smooth tube in configuration CO 
(which is identical to the configuration used by Perry 
[Ref. 2]). In Figure 5.1, the three runs completed are in 
very good agreement with each other. At lower heat fluxes 
( < 4 kW/m^) some scatter is present which disappears at mid- 
range heat fluxes (4-8 kW/m^) . Perry [Ref. 2] found the same 
phenomenon and attributed the scatter to convection regime 
thermal gradients where mixing is inefficient within the 
liquid pool, after discovering large local variations in the 
tube wall and pool thermocouples of up to 5%. In this study, 
local variations in the tube wall and pool thermocouples were 
found to be up to 7% at lower heat fluxes and decreasing as 
heat flux increased, which supports Perry's theory. In 
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Heat Flux/(W/m'2) 




Wall Superheat/(K) 



Figure 5.1 Increasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration CO 
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natural convection, the three data runs started to the left of 
the Churchill and Chu correlation’ [Ref. 9] (shown as a solid 
line) at low heat fluxes and converged toward the Churchill- 
Chu line as the heat flux was increased. This slight shift to 
the left, although consistent during the three runs, was 
within the uncertainty in the region. Good repeatability was 
also found at the onset of nucleate boiling (indicated by 
arrows) and in the nucleate boiling region. Comparing the 
repeatability runs to the representative data set of Perry 
[Ref. 2], at low heat fluxes the data vary slightly but are 
within the uncertainty band. The data merge up to the point 
of incipience, but why Perry's run nucleated earlier than the 
present data is unclear since the runs were, in essence, 
identical. In the boiling region. Perry's data is slightly to 
the left of the present data but again lie within the 
uncertainty band for that region. 

In Figure 5.2, good repeatability is again shown for all 
heat fluxes. Comparing the three present repeatability runs 
to the representative data set of Perry [Ref. 2], again shows 
good reproducibility throughout. For heat fluxes < 1 IcW/m^, 
the slopes of the data curves change, which indicates the 
transition from nucleate boiling bade to natural convection. 
Based on Figures 5.1 and 5.2, data reproducibility is 
considered satisfactory with the single smooth tube worlc 
carried out by Perry [Ref. 2] on the same apparatus. 
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Figure 5.2 Decreasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration CO 
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Figure 5.3 is a comparison of decreasing heat flux 
repeatability runs using an enhanced tube (TURBO-B) in 
configuration CO with the single tube data of Bertsch 
[Ref. 1]. For the present data, good repeatability is again 
obtained throughout. In comparing the current data with the 
data collected by Bertsch [Ref. 1], different boiling curves 
are easily seen. This was a result of two different TURBO-B 
test tubes being used. After examining both test tubes under 
a microscope (using 50x and lOOx magnification) , the TURBO-B 
test tube used in this study had wider notches between the 
'mushroom-like pedestals' (see Figure 3.4), indicating that a 
different fabrication process was used to make the tube. 
Above 20 kW/m^ on the plot, the current data curves are to the 
left of Bertsch' s data (about 2 °C less wall superheat), 
indicating better heat transfer. It is believed this is due 
to the wider notches in the tube surface allowing more bubbles 
to depart the tube surface at higher heat fluxes, removing 
more heat than the tube used by Bertsch. At about 10 kW/m^, 
the trend is reversed, the current data crosses and continues 
to fall below Bertsch 's data (about 0.3 °C more wall 

superheat) . At low heat fluxes, the wider notches allow the 
refrigerant to flood back into the boiling channels (see 

Figure 3.4) and reduce the number of active sites, reducing 

the heat transfer from the tube. 

Figure 5.4 is a comparison of increasing heat flux 

repeatability runs using smooth tube configuration Cl with the 
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Figure 5.3 Decreasing Heat Flux Repeatability Runs using 
TURBO-B Tube Configuration CO 
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Figure 5.4 Increasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration Cl with the Lower Tube Unheated 
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lower tube unheated and a single smooth tube data set 
(D0728ASC0) from Figure 5.1. The plot shows some scatter of 
data at low heat fluxes that is attributed to the liquid pool 
temperature gradients, as discussed previously. Again, good 
repeatability is found in the mid-range heat flux and boiling 
regions. In comparison to the single tube data (configuration 
CO) and the Churchill and Chu correlation [Ref. 9], the upper 
tube heat transfer performance of the Cl configuration is 
clearly enhanced in natural convection. This enhancement is 
probably due to the liquid pool surface being much closer to 
the upper test tube (approximately 10 mm) when it is position 
Cl, increasing the circulation (i.e., local fluid velocity) 
around the upper tube and ultimately increasing heat transfer 
performance of the tube. Lake [Ref. 3], in his comparison of 
single tube R-114 data to his Cl configuration data, found the 
same enhancement. He concluded that the lower tube must be 
modifying the flow over the upper tube. In actuality, it is 
probably a combination of both the decreased distance to the 
pool surface and the lower tube modifying the flow over the 
upper tube that causes the enhancement. From this plot, it 
was concluded that the program DRPJY was performing properly. 

Figure 5.5 is the corresponding decreasing heat flux data 
runs for Figure 5.4. Good repeatability is again found for 
the Cl configuration runs. Comparing the single tube data 
(D0728ASC0 from Figure 5.2) to the Cl configuration data. 
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Figure 5.5 Decreasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration Cl with the Lower Tube Unheated 
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indicates that in nucleate boiling with decreasing heat 
fluxes, both tube configurations perform equally well. 

Figure 5.6 is the decreasing heat flux repeatability runs 
using TURBO-B tubes in configuration Cl with the lower tube 
unheated. Good repeatability is evident. Comparing the 
single TURBO-B data set (D1012TBC0) from Figure 5.3 to the Cl 
data, at high heat fluxes both tube configurations perform 
equally well, as was for the smooth tube case (Figure 5.5). 
At low heat fluxes (< 8 )cW/m^) , the Cl configuration data is 
to the left of the CO configuration data, indicating an 
enhancement due to the addition of the lower tube and moving 
the test tube closer to the surface of the liquid pool; this 
was not seen for the smooth tube case (Figure 5.5). As will 
be discussed later, the addition of a lower tube heat flux, 
shifts the Cl data curve even more towards the left. 

Based on Figures 5.1 through 5.6, the repeatability of the 
current data and reproducibility of previous single tube worlc 
performed by Perry [Ref. 2], is considered satisfactory. 

B. INFLUENCE OF THE LOWER TUBE 

To investigate the effects of a lower heated tube on an 
upper tube, five different tube configurations and fixed lower 
tube heat flux settings were used. The five tube 
configurations selected were P/D ratios of 2.0, 1.8, 1.5, 1.2 
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Figure 5 . 6 Decreasing Heat Flux Repeatability Runs using 
TURBO-B Tube Configuration Cl with the Lower Tube Unheated 
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and 1.2 (rotated 30°) . P/D of 2.0 was selected because it 
matches not only the vertical pitch used in naval shipboard 
flooded evaporators, but also the pitch u..ed in the flooded 
evaporator bundle apparatus at the Naval Postgraduate School 
Heat Transfer Laboratory. The fixed lower tube heat flux 
settings (LTHFS) used were 0, 1, 5, 10 and 20 kW/m^ . In 
addition to these runs, data runs for each tube configuration 
were also done using a lower tube heat flux setting that was 
equal to the upper tube heat flux setting (UTHFS) . Finally, 
at a fixed UTHFS of 5, 10 and 20 kW/m^, the LTHFS was varied 
to investigate the influence of the strength of the convection 
effects from below. For each data run, the same tube surface 
was used for the upper and lower tube; mixing of the tube 
types was not allowed . For the smooth tube runs, the upper 
tube was started in the natural convection region at low heat 
fluxes, was run up and down through various heat flux 
settings, and was returned to the natural convection region. 
For the TURBO-B tubes, only decreasing upper tube heat flux 
runs were recorded to prevent overlapping/repeating data (the 
TURBO-B tubes did not enter the natural convection region at 
low heat fluxes) . 

To evaluate the effect of a lower heated tube on an upper 
tube, the data is presented by the effect of changing the P/D 
ratio for a common LTHFS, the effect of varying the LTHFS for 
a common P/D ratio and the effect of varying the lower tube 
heat flux for a fixed upper tube heat flux setting. 
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1 . The Effect of Tube Spacing 

a. Smooth Tabes vlth Lover Tube Unheated 

Figure 5.7 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with the lower 
tube unheated and the corresponding Churchill and Chu 
correlation [Ref. 9]. As seen in the repeatability runs, at 
low heat fluxes the data is shifted to the left of the 
Churchill-Chu line and as the heat flux is increased, the data 
tends to approach the Churchill-Chu line up to the incipient 
boiling point (indicated by arrows) . This implies that the 
circulation occurring within the evaporator enhances the 
natural convection that occurs from a single tube immersed in 
an infinite expanse of fluid. Following nucleation of the 
upper tube, the boiling curve for P/D of 1.8 is to the left of 
the others and appears to be yielding the best heat transfer 
performance. It is unclear why this is so, because it is so 
repeatable and the same smooth test tubes were used. Lake 
[Ref. 3] also found that P/D of 1.8 gave the best heat 
transfer performance. 

Figure 5.8 shows the corresponding decreasing heat 
flux data for Figure 5.7. With decreasing heat flux, a P/D of 
1.8 again provides the best heat transfer performance until 
the lowest heat flux region, as did Lake's data [Ref. 3]. 
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Figure 5.7 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with the Lower Tube Unheated 
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Figure 5 . 8 Conparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with the Lower Tube Unheated. 
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Jb. Smooth Tubes with LTHFS of 1 kW/m^ 

Figure 5.9 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with a heat 
flux of 1 kW/m^ on the lower tube. Both the upper and lower 
tubes for all configurations started in the natural convection 
region. At low heat fluxes, some scatter is present which 
becomes smaller as the data approaches the incipient boiling 
points and tends towards the Churchill and Chu correlation 
[Ref. 9], as seen previously. Note that the data for a P/D of 
1.2 crosses the Churchill-Chu line and appears to nucleate 
earlier than the others. For this small spacing, the heated 
plume from the lower tube can actually reduce natural 
convection from the top tube since the local ambient pool 
temperature is increased. This creates a larger wall 
temperature and a larger wall superheat (Twall - Tsat) . After 
nucleation, a P/D of 1.8 again provides better heat transfer 
in the boiling region than the other P/D ratios. This trend 
suggests that in the boiling region the enhancement provided 
by a P/D of 1.8 is due to the tube spacing itself and not from 
heating the lower tube. 

Figure 5.10 is the decreasing heat flux data for to 
Figure 5.9. A P/D of 1.8 provides the best heat transfer 
performance until approximately 4 kW/m^, which is where the 
boiling curves seem to merge. With a LTHFS of 1 kW/m^, it 
appears to force the decreasing boiling curves to merge at a 
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Figure 5 . 9 Coirparisoa of Smooth Tube P/D Ratios for 
Increasing Heat Flux with 1 kW/m^ on the Lower Tube 
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Figure 5.10 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 1 kW/m^ on the Lower Tvibe 
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higher heat flux (4 kW/m^) than with the lower tube unheated 
(< 1 kW/m^) . 

c. Smooth Tubes vith LTHFS of 5 kW/nt^ 

Figure 5.11 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with 5 kW/m^ on 
the lower tube. Both upper and lower tubes for all the 
configurations started in the natural convection region. At 
low heat fluxes, the data is to the left of the Churchill and 
Chu correlation [Ref. 9], as seen in Figures 5.7 and 5.9. As 
the heat flux is increased, some interesting things happen. 
The data for a P/D of 1.2 (rotated 30°) quickly crosses the 
Churchill-Chu line and begins to parallel it, then nucleates 
at a low heat flux of about 3 kW/m^ . This is believed to be 
an anomaly and not an effect. Of the remaining P/D ratios, 
1.8 is the only one to tend toward the Churchill-Chu line with 
increasing heat flux. The P/D ratios of 2.0, 1.5 and 1.2 seem 
to parallel the Churchill-Chu line, probably due to the LTHFS 
intensifying the forced convection around the upper tube. In 
all cases, the lower tube remained in natural convection until 
the upper tube nucleated. Once the upper tube nucleated, the 
lower tube also nucleated. In the boiling region, a P/D of 
1.8 again gave the best enhancement. 

Figure 5.12 is the decreasing heat flux data for 
Figure 5.11. Above 10 kW/m^, P/D of 1.8 is clearly providing 
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Figure 5.11 Coitiparison of Smooth Txibe P/D Ratios for 
Increasing Heat Flux with 5 kW/m^ on the Lower Tube 
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Figure 5.12 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 5 kW/m^ on the Lower Tube 
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the best heat transfer. Below 10 kW/m^, the upper tube is 
returning to natural convection, but the lower tube is still 
nucleating. The bubbles departing the lower tube impinge upon 
the upper tube causing a significant increase in heat transfer 
performance, from bubble sweeping, compared to the data for a 
LTHFS of 0 and 1 kW/m^ Also, a LTHFS of 5 kW/m^ eliminates 
any upper tube hysteresis effect. 

d. Smooth Tubes v±th LTHFS of 10 kW/a^ 

Figure 5.13 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with 10 kW/m^ 
on the lower tube. For this LTHFS, the lower tube nucleated 
at the start of the run for a P/D of 1.5, but the upper tube 
remained in natural convection. The effect of the bubbles 
from the lower tube impacting and sweeping over the upper tube 
is seen to significantly enhance the heat transfer 
performance of the upper tube in the low heat flux regions. 
For the other P/D ratios, the incipient boiling points were 
very random; at the largest spacing, nucleation occurred 
quickly, whereas at the smallest spacing, nucleation was 
delayed. Possibly, the largest spacing (P/D of 2.0) allowed 
the heated plume from the lower tube to thin the boundary 
layer around the upper tube, but the closest spacing (P/D of 
1.2) enabled the heated plume to 'skirt' around the upper tube 
and not disturb the upper tube boundary layer as much, 
delaying the point of incipience. Again, once the upper tube 
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Figure 5.13 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with 10 kW/m^ on the Lower Tube 
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nucleated, the lower tube also nucleated shortly thereafter. 
In the fully developed nucleate boiling region, a P/D of 1.8 
still gives the best performance. 

Figure 5.14 is the decreasing heat flux data for 
Figure 5.13. A P/D of 1.8 provided the best heat transfer 
until the heat flux was lowered to about 15 kW/m^, where no 
gain over the other P/D ratios is noticed. This concurs with 
the findings of Lake [Ref. 3] . Below 10 kW/m^, the nucleating 
lower tube is definitely enhancing the upper tube performance 
for all P/D ratios, although, the P/D of 1.2 (rotated 30°) is 
slightly to the right of the other curves. This is probably 
caused by the angular offset of the upper tube, not allowing 
the bubbles from the lower tube to directly impinge upon the 
upper tube. Therefore, a 30° rotation slightly decreases the 
performance of the upper tube compared to an in-line tube 
configuration . 

e. Smooth Tubes with LTHFS of 20 kW/a^ 

Figure 5.15 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with 20 kW/m^ 
on the lower tube. In each run, the lower tube was nucleating 
at the beginning of the run and remained nucleating 
throughout. The only effect from P/D spacing is noticed above 
20 kW/m^ where a P/D of 1.8 again provides the best 
performance . 
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Figure 5.14 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 10 kW/m^ on the Lower Tube 
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Figure 5.15 Con^arison of Smooth Ttibe P/D Ratios for 
Increasing Heat Flux with 20 kW/m^ on the Lower Tube 
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Figure 5.16 is the decreasing heat flux data for 
Figure 5.15. Again, below 20 kW/m^, only a P/D of 1.2 
(rotated 30°) drops slightly below the other curves at mid- 
range heat fluxes. It appears that this LTHFS enables more 
bubbles from the lower tube to impinge upon the upper tube, 
compared to a LTHFS of 10 kW/m^. 

f. Smooth Tx2bes vlth UTHFS equal LTHFS 

Figure 5.17 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with the UTHFS 
equal to the LTHFS throughout the entire length of the run. 
At the beginning of all the runs, both tubes were in the 
natural convection region. In the natural convection region, 
a P/D of 1.2 (rotated 30°) is the only configuration that is 
very close to the Churchill and Chu correlation [Ref. 9] . 
This was not expected nor can be explained. Again, a random 
incipient point was obtained and in the boiling region, a P/D 
of 1.8 was the best performer. 

Figure 5.18 is the decreasing heat flux data for 
Figure 5.17. Above 20 kW/m^, a P/D of 1.8 still provides the 
best performance and below 20 kW/m^ it appears that as the P/D 
decreases, so does the enhancement of the upper tube. Between 
9 and 5 kW/m^, it's believed that the tubes shift from the 
boiling region back to natural convection. This may cause the 
decrease of enhancement because, as speculated earlier, as the 
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Figure 5.16 Con^sarison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 20 kW/m*” on the Lower Tube 
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Figure 5.17 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with UTHFS equal to LTHFS 
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Figure 5.18 Coitparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with UTHFS equal to LTHFS 
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P/D is decreased so too are the convection effects from the 
lower tube. 

g. TURBO-B Tubes vith Lover Tube Unheated 

Figure 5.19 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with lower tube 
unheated. At the start of each run, an undetermined amount of 
nucleation was occurring from the lower tube due to the 
vigorous boiling action off the upper tube at the run start 
heat flux (100 kW/m^) . A very consistent boiling curve above 
10 kW/m^ indicates that there is no enhancement from the lower 
tube at any P/D ratio. Below 5 kW/m^, the rotated 

configuration data agrees well with the single tube data in 
Figure 5.3. This indicates that the lower tube is not 

modifying the flow around the upper tube and the upper tube is 
performing as a single tube which provides less heat transfer 
than the in-line configurations. 

h. TUBBO-B Tubes vith LTHFS of 1 kW/m^ 

Figure 5.20 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with 1 kW/m^ on 
the lower tube. At the start of each run, both the upper and 
lower tubes were nucleating. Very little change is noticed 
between Figures 5.19 and 5.20. Above 20 kW/m^, the plots are 
identical. Below 10 kW/m^, the LTHFS of 1 kW/m^ moved the data 
curves slightly to the left. This indicates that the greater 
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Figure 5.19 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with the Lower Tube Unheated 
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Figure 5 . 20 Con^arison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with 1 kW/m^ on the Lower Tube 
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the LTHFS, the more enhanced the mid-range and lower heat flux 
region becomes. As seen in the previous figure, a P/D of 1.2 
(rotated 30°) is again less enhanced than the other P/D 

ratios . 

i. TUPBO-B Tubes with LTHFS of 5 JcW/m^ 

Figure 5.21 shows the decreasing heat flux data 

talcen for the five tube configurations (C1-C5) with 5 )cW/m^ on 

the lower tube. At the start of each run, both tubes were 
vigorously boiling. Above 20 )cW/m^, the data is consistent 
with the two previous plots, but below 10 kW/m^, a definite 
separation between a P/D of 1.2 (rotated 30°) and the other 
P/D ratios is seen. Also for heat fluxes < 10 kW/m^, the in- 
line configurations provide better heat transfer than the same 
configurations seen in Figure 5.20. This indicates that as 
LTHFS increases, more bubbles depart the lower tube and 

increase the upper tube heat transfer performance for the in- 
line configurations. For the rotated configuration, the 
bubbles must be departing the lower TURBO-B tube differently 
(compare to a lower smooth tube) and are not impinging upon 
the upper TURBO-B tube. This explains the approximately 40 - 
45% decrease in heat transfer performance for the rotated 
configuration compared to the in-line configurations. 

j. TUPBO-B Tubes with LTHFS of 10 kW/m^ 

Figure 5.22 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with 10 kW/m^ 
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Figure 5.21 Conparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with 5 kW/m^ on the Lower Tube 
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on the lower tube. Again, both tubes were nucleating at the 
beginning of each run. There are no obvious changes above 20 
kW/m^ from the previous plots. Below 10 kW/m^, all the curves 
seem to have shifted to the left compared to the previous 
case . 

k. TUBBO-B Tubes vith LTHFS of 20 kW/nf 

Figure 5.23 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with 20 kW/m^ 
on the lower tube. Figures 5.22 and 5.23 show similar 
results, therefore, increasing the LTHFS from 10 to 20 kW/m^ 
provides a small change in upper tube heat transfer 
performance . 

l. TUBBO-B Tubes vith UTHFS equal LTHFS 

Figure 5.24 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with the UTHFS 
equal to the LTHFS. There is not much change in this plot 
relative to the trend set by the previous plots. 

2. The Effect of a Fixed Lower Tube Heat Flux 
a. Smooth Tubes vith P/D of 2.0 

Figure 5.25 shows the increasing heat flux data 
taken for a P/D of 2.0 at five different heat flux settings on 
the lower tube. The Churchill and Chu correlation [Ref. 9] is 
also shown for comparison. For a LTHFS of 20 kW/m^ the lower 
tube had nucleated immediately at the beginning of the run and 
remained nucleating for the length of the run. All other 
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Figure 5.23 Conparison of TURBO-B Tube P/D Ratios for 
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Figure 5 . 25 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 2.0 
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LTHFS data runs, started in natural convection. In the 
natural convection region, the data are shifted to the left of 
the Churchill and Chu correlation [Ref. 9], but approach it as 
heat flux is increased. This is probably caused by strong 
circulation patterns within the small evaporator enhancing the 
heat transfer mentioned previously. In natural convection, 
the effect of a LTHFS had little influence on the upper tube 
heat transfer behavior; in fact, this is repeated for all the 
P/D ratios. The point of incipience was very repeatable, with 
the exception of a LTHFS of 10 kW/m^. A heat flux of 10 kW/m^ 
is approximately the heat flux that 'trips' the nucleation 
sites causing boiling to occur from a smooth tube, so shortly 
after the run started, the lower tube nucleated and caused the 
upper tube to nucleate quickly. At high heat fluxes, all the 
data falls on one nucleate boiling curve indicating that the 
LTHFS has little effect on the heat transfer from the upper 
tube in this region. 

Figure 5.26 shows the decreasing heat flux data for 
the data run for Figure 5.25. Important to note that the 
upper tube starts and finishes each run in natural convection, 
but the lower tube might not. For a LTHFS of 0 and 1 kW/m^, 
the lower tube starts and finishes in natural convection; for 
a LTHFS of 20 kW/m^, it starts and finishes in nucleate 
boiling; and for a LTHFS of 5 and 10 kW/mS it starts in 
natural convection and finishes in nucleate boiling. At high 
heat fluxes, all the data lay on one nucleate boiling curve as 
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Figure 5.26 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 2.0 



91 



expected. Below 20 kW/m^, five separate curves are seen. The 
curves with a LTHFS of 0 and 1 kW/m^ return to the natural 
convection region, indicating a deactivation of the nucleation 
sites as the upper tube heat flux is decreased, showing the 
hysteresis effect of the upper tube. At a heat flux below 5 
kW/m^, the curves with a LTHFS of 5, 10 and 20 kW/m^ begin to 
parallel the natural convection curve found on the increasing 
heat flux plot and have significant enhancement over the 0 and 
1 kW/m^ curves at low heat fluxes. Also, the hysteresis 
effect is eliminated by the vigorously nucleating lower tube. 
This same phenomena was found by Lake [Ref. 3] . He attributed 
this enhancement to strong convection effects on the upper 
tube from the bubbles nucleating off the lower tube. 
Therefore, significant enhancement of heat transfer of an 
upper tube is possible when the lower tube is nucleating and 
very little enhancement is noticed when both tubes are in 
natural convection. 

Jb. Smooth Tubes with P/D of 1.8 

Figure 5.27 shows the increasing heat flux data 
taken for a P/D of 1.8 at five lower tube heat flux settings. 



For a LTHFS of 


20 


kW/m^, the 


lower tube 


had 


nucleated 


immediately at 


the 


beginning 


of the run 


and 


remained 



nucleating for the entire length of the run. All other LTHFS 
started in natural convection. In the natural convection 
region, the data are shifted to the left of the Churchill and 
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Figure 5.27 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 1.8 
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Chu correlation [Ref. 9], but approach it as heat flux is 
increased, as was the case with a P/D of 2.0. Again, the 
effect of the LTHFS had little influence of the upper tube 
heat transfer behavior. The point of incipience was again, 
very repeatable. This time, unlike for a P/D of 2.0, the 
lower tube remained in natural convection until the upper tube 
nucleated. Again, at high heat fluxes, all the data falls on 
one nucleate boiling curve indicating that the LTHFS has 
little effect upon the heat transfer in this region. 

Figure 5.28 shows the decreasing heat flux data for 
the data runs for Figure 5.27. The data trend is similar to 
that found in Figure 5.26. Therefore, significant enhancement 
of heat transfer from an upper tube is again seen when the 
lower tube is nucleating with larger than 5 kW/m^ on the lower 
tube . 

c. Smooth Tubes with P/D of 1.5 

Figure 5.29 shows the increasing heat flux data 
taken for a P/D of 1.5 at five lower tube heat flux settings. 
For a LTHFS of 10 and 20 kW/m^, the lower tube had nucleated 
immediately at the beginning of the run and remained 
nucleating for the length of the run. All other LTHFS started 
in natural convection. This plot shows similar results seen 
on Figure 5.27 and most of the comments still apply. 
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Figure 5 . 28 Conparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Txibes with P/D of 1.8 
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Figure 5.29 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Txibes with P/D of 1.5 
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Figure 5.30 shows the decreasing heat flux data for 
the data runs shown in Figure 5.29. Again, similar results 
are obtained as before. Significant enhancement in the heat 
transfer from an upper tube is obtained when the lower tube is 
nucleating with a heat flux larger than 5 kW/m^. 

d. Smooth Tubes with P/D of 1.2 

Figure 5.31 shows the increasing heat flux data 
taken for a P/D of 1.2 at five lower tube heat flux settings. 
For a LTHFS of 20 kW/m^, the lower tube had nucleated 
immediately at the beginning of the run and remained 
nucleating for the length of the run. All other LTHFS started 
in natural convection. This plot shows similar results seen 
on Figure 5.29 and most of the comments still apply. 

Figure 5.32 shows the decreasing heat flux data for 
the data runs shown in Figure 5.31. This plot shows similar 
results seen on Figure 5.30 and again, most of the comments 
still apply. 

e. Smooth Tubes vith P/D of 1.2 (Rotated 30°) 

Figure 5.33 shows the increasing heat flux data 

taken for a P/D of 1.2 (rotated 30°), at five lower tube heat 
flux settings. For a LTHFS of 20 kW/m^, the lower tube had 
nucleated immediately at the beginning of the run and remained 
nucleating for the length of the run. All other LTHFS started 
in natural convection. The increasing data for this 
configuration is more concentrated than the others seen 
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Figure 5.30 Conparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Snusoth Tubes with P/D of 1.5 
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Figure 5.31 Conparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 1.2 
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Figurs 5.33 Coitpsifison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 1.2 
(Rotated 30'*) 
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previously. This might be due to the offset orientation of 
the upper tube causing the data to be less scattered. 

Figure 5.34 shows the decreasing heat flux data for 
the data runs shown in Figure 5.33. At high heat fluxes, all 
the data lay on one nucleate boiling curve. Below 10 kW/m^, 
separate curves are formed and as seen before. Significant 
enhancement of heat transfer from an upper tube is obtained 
when the lower tube is nucleating with a heat flux larger than 
5 kW/m^ 

t. TUBBO-B Tubes with P/D of 2.0 

Figure 5.35 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
At the start of each run, a very consistent boiling curve is 
seen above 20 kW/m^, indicating that above a UTHFS of 20 kW/m^, 
there is no enhancement from increasing the LTHFS. Below 10 
kW/m^, raising the LTHFS does increase the heat transfer 
performance of the upper tube. At an upper tube heat flux of 
2 kW/m^, an increase in performance of 2.4 is found going from 
an unheated lower tube to a LTHFS of 20 kW/m^, compared to an 
increase of 3.6 for the smooth tube array (Figure 5.26). 
gr. TURBO-B Tubes with P/D of 1.8 

Figure 5.36 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
Similar results are found compared to the P/D of 2.0 case. At 
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Figure 5.34 Coniparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 1.2 
(Rotated 30 *) 



103 



Upper Tube Heat Flux/(W/m‘2) 



10 ' 












A 



s 






q^Ox- 



10 ^ _ 



<sP >t 



03 



0 X 



^ 0 



x+ 



10^ ^ 



+ q" lower tube - 0 kW/m‘2 

X q* lower tube - 1 kW/m‘2 

0 q* lower tube - 5 kW/m*2 

° q* lower tube - 10 kW/m‘2 

A q" lower tube - 20 kW/m‘2 



10 



-1 



10 ° 10 ^ 
Wall Superheat/(K) 



Figure 5.35 Coirparison of Lower Txibe Heat Flux Settings on 
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Figure 5.36 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.8 
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an upper tube heat flux of 2 kW/m^, an increase in performance 
of 1.6 is found going from a unheated lower tube to a LTHFS of 
20 kW/mS compared to an increase of 3.7 for the smooth tube 
array (Figure 5.28) . Unlike the smooth tube case where a P/D 
of 1.8 produced the best heat transfer performance in nucleate 
boiling, no optimal P/D ratio was found for the TURBO-B tube 
array. 

h. TUPBO-B Tubes vith P/D of 1.5 

Figure 5.37 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
Similar results are found compared to the P/D of 1.8 case. At 
an upper tube heat flux of 2 kW/m^, an increase in performance 
of 2.2 is found going from a unheated lower tube to a LTHFS of 
20 kW/m^, compared to an increase of 3.5 for the smooth tube 
array (Figure 5.30). 

i. TURBO-B Tubes vith P/D of 1.2 

Figure 5.38 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
Similar results are found compared to the P/D of 1.5 case. At 
an upper tube heat flux of 2 kW/m^, an increase in performance 
of 2.2 is found going from a unheated lower tube to a LTHFS of 
20 kW/m^, compared to an increase of 5.4 for the smooth tube 
array (Figure 5.32). 
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Figure 5.37 Comparison of Lower T\ibe Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.5 
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Figure 5.38 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.2 
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j. TURBO-B Tubes with P/D of 1.2 (Rotated 30°) 

Figure 5.39 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
A very consistent boiling curve is again seen above 10 kW/m^. 
Below 10 kW/m% the data seems less scattered and slightly to 
the left of any other configuration, indicating that rotating 
the upper tube reduces the enhancement produced from 
increasing the LTHFS in the low heat flux regions. At an 
upper tube heat flux of 2 kW/m^, an increase in performance of 
1.7 is found going from a unheated lower tube to a LTHFS of 20 
kW/m^, compared to an increase of 3.5 for the smooth tube 
array (Figure 5.34). 

3 . The Effect of a Varieible Lower Tube Heat Flux 
a. Smooth Tubes vith P/D of 2.0 

Figure 5.40 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is varied. The three fixed UTHFS used were 5, 10 
and 20 kW/m^. The labeling in the plot legend describes which 
data were taken using increasing or decreasing lower tube heat 
flux. The sharp transitions in the increasing curves are the 
points where the lower tube nucleated and caused the upper 
tube to shift from natural convection to nucleate boiling 
(visually verified during the data run) , significantly 
increasing the upper tube heat transfer coefficient, as 
expected. The maximum upper tube heat transfer coefficient 
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Figure 5.39 Comparison of Lower Tvibe Heat Flux Setting on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.2 
(Rotated 30®) 
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Figure 5 . 40 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 2 . 0 
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for all three UTHFS was 1850 W/m^-K, which indicates that at 
a lower tube heat flux of 80 kW/m^, the heat transfer 
performance of the upper tube is independent of the lower tube 
heat flux. The decreasing curves show that for an UTHFS of 20 
kW/m^, the lower tube heat flux has little influence on the 
upper tube heat transfer. For the UTHFS of 10 kW/m^, as the 
lower tube heat flux is decreased, the upper tube remains in 
nucleate boiling, but for the UTHFS of 5 kW/m^, the upper tube 
returns to the natural convection region as the lower tube 
heat flux is decreased. 

jb. Smooth Tubes vith P/D of 1.8 

Figure 5.41 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is varied. Similar results were obtained as in the 
P/D of 2.0 case, except that the curves did not merge at the 
highest heat flux setting. The maximum heat transfer 
coefficient for the UTHFS of 5 and 10 kW/m^ converged at 
approximately 1700 W/m^-K and for the UTHFS of 20 kW/m^ the 
value was about 1900 W/m^-K. It is not known why the curves 
did not merge at high lower tube heat fluxed. Comparing these 
results to the P/D of 2.0 results, at the smallest lower tube 
heat fluxes, the P/D of 1.8 provides better heat transfer. 
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Figure 5 . 41 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.8 
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Figure 5 . 42 Heat TrcUisfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.5 
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c. Smooth Tubes vith P/D of 1.5 

Figure 5.42 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is varied. The maximum heat transfer coefficient 
for all three of the UTHFS converged at approximately 1800 
W/m^*K, similar to the P/D of 2.0 case. The decreasing curves 
show very similar behavior found at a P/D of 2.0. Comparing 
these results to data of a P/D of 2.0 and 1.8, a P/D of 2.0 
and 1.5 are similar in the amount of enhancement achieved, 
whereas a P/D of 1.8 outperforms them both, in the mid to low 
heat flux regions. 

d. Smooth Tubes vith P/D of 1.2 

Figure 5.43 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is increased. The maximum heat transfer coefficient 
for all three of the UTHFS converged, as did the P/D of 2.0 
and 1.5, at approximately 1990 W/m^’K, highest value of all 
the configurations tested. Therefore, the best heat transfer 
performance was given by a P/D ratio of 1.8 when the lower 
tube heat flux was held constant, but a P/D ratio of 1.2 when 
the upper tube heat flux was held constant. No decreasing 
curves were recorded for these data runs. 

e. Smooth Tubes vith P/D of 1.2 (Rotated 30°) 

Figure 5.44 shows the heat transfer performance of 

an upper tube, held at a constant heat flux, as the lower tube 
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Figure 5 . 43 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 
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Figure 5 . 44 Heat Tremsfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 (Rotated 30*) 
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heat flux is increased. As seen for the P/D of 1.8 case, the 
curves do not merge at the highest lower tube heat flux. The 
maximum heat transfer coefficient for the two lower curves 
converged at about 1500 W/m^-K and the upper curve was 1725 
W/m^‘K. Comparing Figures 5.43 and 5.44, for an UTHFS of 20 
kW/m^, the data curves are identical for lower tube heat 
fluxes up to 8 kW/m^, indicating that the forced convection 
effects from the lower tube have very little influence on the 
upper tube heat transfer behavior. Above lower tube heat 
fluxes of 10 kW/m^, the in-line configuration provides better 
heat transfer performance because the lower tube is in 
nucleate boiling and all of the bubbles are impacting the 
upper tube providing enhancement by bubble pumping; as 
opposed to the offset configuration, where most of the bubbles 
coming off the lower tube miss the upper tube. The same idea 
holds true for the other data curves. The in-line 
configurations are more enhanced due to bubble pumping from 
the lower tube. No decreasing curves were recorded for these 
data runs. 

f. TUPBO-B Tubes vith P/D of 2.0 

Figure 5.45 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. At high lower tube heat fluxes, all 
three of the UTHFS appear to perform almost equally. For an 
UTHFS of 20 kW/m^, the lower tube heat flux has little 
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Figure 5 . 45 Heat Transfer Performance of a TUFBO-B Upper 
Tube at a Constamt Heat Flux for a P/D of 2.0 
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influence on the upper tube heat transfer and for an UTHFS of 
5 and 10 kW/m^, as the lower tube heat flux is decreased so 
follows the performance of the upper tube (more for 5 than 10 
kW/m^) . 

g. TUPBO-B Txibes vith P/D of 1.8 

Figure 5.46 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. Comparing these curves to the ones 
for a P/D of 2.0, a small enhancement in performance is gained 
in the mid and low heat flux regions using a P/D of 1.8, which 
is an expected result. 

h. TUPBO-B Txjbes vith P/D of 1.5 

Figure 5.47 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. There is very little difference in 
performance between a P/D of 1.5 and 1.8. Comparing the P/D 
of 1.5 to the 2.0 data, a small decrease in performance is 
noted for lower tube heat fluxes > 20 kW/m^ using the P/D of 
1.5 configuration, which cannot be explained. 

i. TUPBO-B Tubes vith P/D of 1.2 

Figure 5.48 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. Similar performance is noted for a 
P/D ratio of 1.2 and 1.5, therefore, the comments made 
previously still apply. 
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Figure 5.46 Heat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.8 
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Figure 5 . 47 Heat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.5 
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Figure 5.48 Heat Tramsfer Perfonoance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 
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j. TURBO-B Tubes with P/D of 1.2 (Rotated 30°) 

Figure 5.49 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. Comparing Figures 5.48 and 5.49, for 
an UTHFS of 20 kW/mS the data curves are almost identical, 
indicating that the lower tube orientation and heat flux have 
very little influence on the upper tube heat transfer 
behavior. For the other two data curves, the in-line 
configuration provides better heat transfer performance 
because of bubble pumping^ discussed earlier. For the offset 
configuration, all of the data curves do not merge and are 
almost horizontal lines, indicating there is no enhancement by 
bubble pumping. It is believed that most of the bubbles 
coming off the lower tube miss the upper tube because the 
bubbles depart the TURBO-B tube surface at the upper most 
section of the tube, which is different than how bubbles 
depart the smooth tube surface. This phenomena needs to be 
further investigated. No decreasing curves were recorded for 
these data runs. 

C. R-124/R-114 PERFORMANCE COMPARISON 

Figure 5.50 compares the decreasing heat flux data for R- 
124 found in this study with the R-114 data of Lake [Ref. 3] 
for smooth tube configurations listed in the plot legend, with 
the lower tube unheated. In the mid to high heat flux 
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Upper Tube Heat Transfer Coefficient/(W/m‘2*K) 




Lower Tube Heat Flux/(W/m*2) 



Figure 5.49 Heat Tr^msfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 (Rotated 30®) 
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regions, the P/D ratios of 1.8 and 1.5 demonstrated improved 
performance with the R-124, but with a P/D of 2.0, little or 
no improvement was found. In the low heat flux region, the 
uncertainty is too large to make any definitive conclusions. 

Figure 5.51 compares the decreasing heat flux data for R- 
124 with the R-114 data of Lake [Ref. 3] for the smooth tube 
configurations listed in the plot legend, with 10 kW/m^ on the 
lower tube. In the high heat flux region, a P/D of 1.8 and 
1.5 again provided better performance using R-124, but for R- 
114, a P/D of 2.0 displayed slightly better performance. In 
the mid-range heat fluxes, R-114 exhibited improved 
performance in each of the three configurations. Again, in 
the low heat flux region, there was too much uncertainty in 
the data to draw any viable conclusions. 
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VI. 



CONCLUSIONS 



The interaction effects of a lower heated tube on pool 
boiling of R-124 from an upper horizontal tube has been 
studied for tube two types: smooth and TURBO-B. Based on the 
data obtained, the following conclusions are made: 



1. In natural convection, the effect of a lower heated tube 
on the heat transfer from an upper tube was small. 

2. For a smooth tube array, a vigorously nucleating lower 
tube eliminated upper tube hysteresis. 

3. For a smooth tube array in nucleate boiling, P/D of 1.8 
gave the best upper tube heat transfer performance. No 
optimal P/D was found for the TURBO-B tube array. 

4. For both arrays with upper tube heat fluxes < 5 kW/m^, 
the enhancing effect of bubbles from a lower tube was 
dramatic. This enhancement increased as lower tube heat 
flux increased. When upper tube heat fluxes were > 20 
kW/m^, all enhancement disappeared. 

5. For a TURBO-B tube array, a 30 degree offset of the upper 
tube reduced the upper tube heat transfer performance 
(compared to the in-line configurations) . This may 
indicate bubbles depart TURBO-B tubes differently than 
smooth tubes. 

6. For a smooth tube array with the lower tube unheated and 
an upper tube heat flux of > 3 kW/m^, the performance of 
R-124 compared to R-114, improved as the P/D ratio was 
reduced. With a nucleating lower tube (at 10 kW/m^) , 
again the performance of R-124 improved as the P/D ratio 
decreased, but only for upper tube heat fluxes of > 40 
kW/m^ . 
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VII. 



RECOMMENDATIONS 



1. Conduct further studies using TURBO-B tubes to verify 
that bubbles do depart the TURBO-B tubes differently than 
smooth tubes. 

2. Conduct similar tests using different enhanced tubes to 
better understand the influence of a lower heated 
enhanced tube on the heat transfer behavior of an upper 
tube . 

3. Modify the existing evaporator to accept a larger number 
of tubes (three and four smooth or enhanced tubes) and 
conduct tests to study the effect of multiple heated 
lower tubes on the heat transfer performance of an upper 
tube. 

4. Conduct similar tests using multiple tube configurations 
with R-124/oil mixtures to determine the effect of oil on 
the heat transfer performance of an upper tube. 
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APPENDIX A. PROGRAM SETUP JY 



Program SETUPJY was used to evaluate the system 
performance prior to commencing a data collection run. This 
program is written in Hewlett-Packard Basic 5.0 language and 
enables the operator to: 



1. Monitor the coolant sump temperature. 

2. Monitor the evaporator liquid pool average temperature. 

3. Monitor all the thermocouple channel output temperatures. 

4. Monitor the voltage, current and power supplied to the 
upper tube heater (main heater) as well as the heater 
element electrical resistance. 

5. Monitor the voltage, current and power supplied to the 
lower tube heater (auxiliary heater) as well as the 
heater element electrical resistance. 
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lot PROGRAM: SETUP JY 

20 i DATS; MAY 15, 1993 

30! PROGRAMMER: LT DEAN SUGIYAMA 

40! MODIFIED 3Y I*ANNIS LAKE JAN 22, 1992 

50! MODIFIED BY GEORGS D- PERRY 10 APR 1993 

60! MODIFIED BY JOE YUSICIAN MAY 15, 1993 

70 COM /Cc/ C(7) 

80 DATA 0.10086091, 25727. 94369, -767345. 3295, 73025595. 31,. 9247486589, 6. 97688B^* 

11,-2.66192E^13 

90 DATA 3.94078E+14 

100 READ C(») 

110 ON KEY 1,15 GOTO 150 

120 PRINTER IS 1 

150 PRINT USING ”2X, ’’’’SELECT OPTION’””’ 

160 PRINT USING "SX, ””0 -MONITOR SUMP’””’ 

170 PRINT USING ”SX, 1 -MONITOR LIQUID””” 

180 PRINT USING ”6X, ””2-CHBCK THERMOCOUPLES ” ” ” 

190 PRINT USING ”6X, ”” 3-CHECR MAIN HEATER””” 

200 PRINT USING ”6X, ” ”4-CHECK AUX HEATERS””” 

210 PRINT USING ”6X,”” 5-EXIT PROGRAM””” 

230 BEEP 
240 INPUT Ido 

250 IF Ido>5 THEN Ido-5 

260 IF Ido-0 THEN 340 

270 IF Ido-1 THEN 520 

280 IF Ido-2 THEN 690 

290 IF Ido-3 THEN 870 

300 IF Ido-4 THEN 870 

310 IF Ido-5 THEN 1170 

340 PRINT 

350 PRINT ”3UMP TEMPERATURE DBG C ” 

360 PRINT 

370 OUTPUT 709; ”AR AF19 AL19 VR5” 

380 OUTPUT 709; ”AS SA” 

390 3ujn-0 

400 FOR J-1 TO 5 

410 ENTER 709;E 

420 Sum-SuiD't’B 

430 NEXT J 

440 Bave«S\un/5 

450 T9mp*FNTv3V(Bave) 

460 PRINT USING ”4X,MD0.DD”;Tamp 

470 BEEP 

480 PRINT 

490 WAIT 5 

500 GOTO 340 

510! 

520 PRINT 

530 PRINT ”LIQUID TEMPERATURE DBG C' 

540 PRINT 

550 OUTPUT 709; ”AR AF16 AL13 VR5” 

560 3up0 

570 FOR I-l TO 2 

580 OUTPUT 709; ”AS SA” 

590 ENTER 709;E 

600 Suzn— Sum'^B 

610 NEXT I 

620 Bavo-Suin/2 

630 Taiop-FNTvsv(Bava) 

640 PRINT USINO ”4X,MDD.DD”;Tamp 
650 BEEP 

660 WAIT 5 

670 GOTO 610 
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6801 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860! 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 



PRINT 

PRINT "TEMPERATURE DEG C" 

OUTPUT 709; "AR AFOO AL19 VR5" 

FOR I=»l TO 20 
OUTPUT 709; "AS SA" 

Suin*0 

FOR J»1 TO 5 
ENTER 709;S 
3uni*Sujn-t'B 
NEXT J 
Save^Sum/ 5 
Temp^FNTvsv ( Save ) 

PRINT TAB(3);I;TA3(15);Tamp 

NEXT I 

3EEP 

WAIT 5 

GOTO 690 

PRINT 

OUTPUT 709; "AR AP20 AL22 VR5" 

FOR I=»l TO 3 
OUTPUT 709; "AS SA" 

Suni*0 

FOR J»1 TO 5 
ENTER 709;B 
Suni*Sujn+B 
NEXT J 

IS I«1 THEN Volt»Sum/5 
IF I«2 AND Ido»3 THEN 

PRINT "MAKE SURE VOLTAGE BOX IS SET TO MAIN HEATERS" 

Amp- Sum/ 5 
END IF 

IF 1*3 AND Ido»4 THEN 

PRINT "MAKS SURE VOLTAGE BOX IS SET TO AUX HEATERS" 
Amp*Sum/ 5 
END IF 
NEXT I 

Amp- AB S ( Amp * 1 . 9 1 8 2 ) 

Volt»ABS(Volt:"25) 

P ower- Vo 1 t * Amp 
Resistance- Volt /Amp 
PRINT 

"BgRP 

PRINT "VOLTAGE(V) CURRENT(A) RBSISTENCS( Ohms ) POWKR(W)" 
PRINT 

PRINT USING "1X,5(MDDDD. DO, 4X)";Volt, Amp, Resistance, Power 
WAIT 5 
GOTO 870 

RgRP 

PRINT 

PRINT "SETUP COMPLETE" 

END 

DBF FNTvsv(V) 

COM /Cc/ C(7) 

T-C(O) 

FOR r-I TO 7 
T»T*^C(I)*V^I 
NEXT I 

T-T>8 . 626897E-2+T*( 3 • 761199B-3-T"5 • 0689259K-5 ) 

RBTTOIf T 
FIfSND 
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APPENDIX B. PROGRAM DRPJY 



This data acquisition/reduction program is written in 
Hewlett-Packard Basic 5.0 language and is listed on the 
following pages. 
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10 t PIL2 NAME: ORPJY 

20 I DATE: MAY IS, 1993 

30 J REVISED VERSION OF DRP72 FOR R-124 (2 INSTRUMENTED TUBES) 

40 J REVISED BY LT. GEORGE PERRY & LT. JOE YUSICIAN 

so COM /Idp/ Idp 

60 PRINTER IS 1 

70 CALL Select 

80 BEEP 

90 PRINT 

100 PRINT "DATA COLLECTION COMPLETE" 

110 END 

1201 

130 SUB Main 

140 COM /I prop/ Ift 

ISO COM /Idp/ Idp 

160 COM /Cc/ C(7),Ical 

170 COM /Wll/ D2,Di,Do,L,Lu,Kcu 

180 DIM Emf(20),T(20),Dla(13),D2a(13),Dla(13),Doa(13),La(13),Lua(13),Kcua(13), 
Et(19),Tn$(4)[lS) 

190 DATA 0. 10086091, 2S727. 94369, -76734S.829S,7802SS9S. 81 

200 DATA -9247486S89,6.97688B+11,-2.66192E*13,3.94078E+14 

210 READ C(") 

220 DATA Smooth, High Flux, Turho-B, High Flux Mod,Turho-B Mod 

230 READ Tn$(") 

240 PRINTER IS 701 

2S0 PRINT 

260 PRINT " Date ;", DATES (TIMEDATE) 

270 PRINT 

280 PRINT USING "10X,""MOTE; Program name : DRPJY""" 

290 BEEP 

3001 INPUT "ENTER INPUT MODE ( 0=30S4A, 1-FILE) ", Im 

310 Iffl-0 

3201 INPUT "SELECT FLUID ( O-R-114, l-R-124) " , Ift 

330 Ift-1 

3401 INPUT "1 OR 2 TUBE OPERATION (ENTER 1 OR 2)",Hwmntu 

3S0 Hwmntu-2 

3601 INPUT "SELECT BEATING MODE (0-ELEC; 1-WATER) ", Ihm 

370 Ihm-0 

3801 INPUT "ENTER THERMOCOUPLE TYPE ( 0-NEW, 1-OLD) ", leal 

390 Ical-0 

400 INPUT "GIVE A NAME FOR THE RAW DATA FILE" ,D2_flleS 

410 PRINT USING "16X, ""New file name: "", lSA";D2_flleS 

420 Slzel-20 

430 CREATE BDAT D2_flle$,Slzel 

440 ASSIGN 9File2 TO D2_flleS 

4S0 1 

460 1 DUMMY FILS UNTIL Nrua KNOWN 
470 Dl_flleS-"DUMMY" 

480 CREATE BDAT Dl_flle$,Slzel 

490 ASSIGN SFllel TO Dl.flleS 

SOO OUTPUT aFllel; Dates 

SIO BEEP 

S20 INPUT "ENTER NR OP DEFECTIVE UPPER TUBE TCS (0-DEFAULT, 2 MAX. )",Idtc 
S30 IF Idtc-0 THEN 

S40 Ldtcl-0 

S50 Ldtc2-0 

S60 PRINT USING "l6X,""No defective UPPER TUBE TCa exist""" 

S70 END IF 

S80 IF Idtc-1 THEN 

S90 BEEP 

600 INPUT "ENTER DEFECTIVE UPPER TUBE TC LOCATION (TEMP l-8)",Ldtcl 

610 PRINT USING "16X, ""UPPER TUBE TC is defective at LOCATION: "" ,DD";Ldtcl 

620 Ldtc2-0 

630 END IP 

640 IF Idtc-2 TBSV 

6S0 BEEP 

660 INPUT "ENTER DEFECTIVE UPPER TUBE TC LOCATIONS (TEMP 1-8 ) " ,Ldtcl ,Ldtc2 

670 PRINT USING "16X, ""UPPER TUBE TC are defective at LOCATIONS; "" ,DD, 4X, DD";L 

dtcl,Ldtc2 
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600 END IF 

690 IF Idtc>2 THEN 

700 BEEP 

710 PRINTER IS 1 

720 BEEP 

730 PRINT "INVALID ENTRY" 

740 PRINTER IS 701 

750 GOTO 510 
760 END IF 

770 OUTPUT aFilel;Ldtcl,Ldtc2 
780 IF HwTnntu=l THEN GOTO 1190 

790 INPUT "ENTER NUMBER OF DEFECTIVE LOWER TUBE TC3 (O^DEFAULT, 2 MAX)",Aldtc 
800 IF Aidtc=0 THEN 
810 Aldtcl=0 
820 Aldtc2=0 

830 PRINT USING "16X,""No defective LOWER TUBE TCs exist""" 

840 END IF 

850 IF Aidtc=l THEN 

860 BEEP 

870 INPUT "ENTER DEFECTIVE TC LOWER TUBE LOCATIONS (TEMP 9-16 ) " , Aldtcl 

880 PRINT USING "16X, ""LOWER TUBE TC is defective at LOCATION DD" ;Aldtcl 

890 Aldtc2=0 

900 END IF 

910 IF Aidtc=2 THEN 

920 BEEP 

930 INPUT "ENTER DEFECTIVE LOWER TUBE TC LOCATIONS (TEMP 9-16) " , Aldtc 1 , Aldtc2 

940 PRINT USING "16X, ""LOWER TUBE TC are defective at LOCATIONS: "" ,DD, 4X,DD";A 

ldtcl,Aldtc2 

950 END IF 

960 IF Aldtc>2 THEN 

970 BEEP 

980 PRINTER IS 1 

990 BEEP 

1000 PRINT "INVALID ENTRY" 

1010 PRINTER IS 701 
1020 GOTO 790 
1030 END IF 

1040 OUTPUT @Fllel;Aldtcl,Aldtc2 
1050 PRINTER IS 1 
1060 BEEP 



1070 


PRINT 


USING 


"2X, 


""Select tube number""" 


1080 


PRINT 


USING 


"6X, 


""4 


Wialand Hard 8 inch""" 


1090 


PRINT 


USING 


"6X, 


w » 5 


HIGH FLUX 8 inch""" 


11001 


PRINT 


USING 


"6X, 




GEWA-K 40 Flns/ln""" 


11101 


PRINT 


USING 


"6X, 


n n "j 


GEWA-K 26 Flns/ln""" 


11201 


PRINT 


USING 


"6X, 


""8 


GEWA-T 19 Flns/ln""" 


11301 


PRINT 


USING 


"6X, 


w « g 


GEWA-T OR GBWA-TY 26 Flns/ln 


11401 


PRINT 


USING 


"6X, 


""10 


THERMOEXCEL-E""" 


11501 


PRINT 


USING 


"6X, 


""11 


THERMOEXCEL-HE " " " 


1160 


PRINT 


USING 


"6X, 


""12 


TURBO-B""" 


11701 


PRINT 


USING 


"6X, 


""13 


GEWA-K 19 Flns/ln""" 


1180 


INPUT 


Itt 









1190 OUTPUT @Filel;Itt 
1200 PRINTER IS 701 

1210 IF Itt=4 THEN PRINT USING " 16X, ""Tubes : Wieland Hard 8 inch""" 

1220 IF Itt’»12 THEN PRINT USING " 16X, " "Tubes : Turbo-B""" 

1230 IF Itt»5 THEN PRINT USING "16X, ""Tubes: High Flux 8 inch""" 

1240! IF Itt<10 THEN PRINT USING "16X,""Tube Number: "",D";Itt 

1250! IF Itt>9 THEN PRINT USING "16X,""Tube Number: "",DD";Itt 

12601 INPUT "ENTER OUTPUT VERSION ( 0» LONG, 1» SHORT, 2 *NONE) ", lov 
1270 lov^O 

12801 INPUT "SELECT (0*LIQ, 1»VAP,2»(LIQ+VAP) /2 ) ", Ilqv 

1290 Ilqv=«2 

1300! 

1310! DIMENSIONS FOR TESTED TUBBS 
13201 

13301 Dl=Diamoter at thormocouplo positions 

1340 DATA .0111125, .0111125, .0111125, .0129540, .012446, .0129540, .0100965 
1350 DATA .0100965, .01157, .01157, .01157, .01157, .01157, .0100965 
1360 READ Dla(*) 

1370 Dl=Dla(Itt) 
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13801 

13901 D2*Dlajneter of teat section to the baae of fins 

1400 DATA .015875, .015875, .015875, . 015824, . 015875 ,. 015824, .01270 

1410 DATA .0127, .0138, .0138, .0138, .0138, .0138, .0127 

1420 READ D2a(*) 

1430 D2=D2a(Itt) 

14401 

14501 Di^Inside diameter of unenhanced ends 

1460 DATA .0127, .0127, . 0 127 ,. 0132 ,. 0 127 , .0132, . 0111125, . 0111125 
1470 DATA .0118, .0118, .0118, .0118, .0118, .0111125 
1480 READ Dia(*) 

1490 Di=Dia(Itt) 

15001 

13101 Do^Outside diameter of unenhanced ends 

1520 DATA .015875, .015875, .015875, . 015824 ,. 015875 ,. 015824 , . 01270, . 01270 
1530 DATA .01331, .01331, .01331, .01331, .0158, .0127 
1540 READ Doa( *) 

1550 Do*Doa(Itt) 

15601 

15701 L=Length of enhanced surface 

1580 DATA .1016, .1016, . 1016, . 1016, .2032, .2032, .2032, .2032, .2032, .2032, .2032, .20 

32. . 2032. .2032 
1590 READ La(*) 

1600 L=La(Itt) 

16101 

16201 Lu=Length of unenhanced surface at the ends 

1630 DATA .0254, .0254, .0254, .0254, .0762, .0762, .0762, .0762, .0762, .0762, .0762, .07 

62. . 0762. .0762 
1640 READ Lua(^) 

1650 Lu*Lua(Itt) 

16601 

16701 Kcu=Thermal Conductivity of tube 

1680 DATA 398,344,344,45,344,45,344,344,398,398,398, 398,398,344 
1690 READ Kcua(*) 

1700 Kcu=Kcua(Itt) 

17101 

1720 Xn=.8 

1730 Fr=.3 

1740 Cf=3.7037E+10 

1750 A*PI * ( Do^2-Di^2 ) / 4 

1760 P=PI-^Do 

1770 Nsub=0 

1780 Ntc=20 

1790 J=1 

1800 Repeat: 1 

1810 Dtld=»2.22 

1820 Ido=2 

1830 ON KEY 1,15 RECOVER 1800 
1840 PRINTER IS 1 

1850 PRINT USING " 2X SELECT OPTION"’*" 

1860 PRINT USING "6X,""0=TAKE DATA""" 

1870 PRINT USING "6X,""1=SET UPPER TUBE HEAT FLUX (MAIN)""" 

1880 PRINT USING "6X,""2»SET Tsat""" 

1890 PRINT USING "6X, ""3-SET LOWER TUBE HEAT FLUX (AUX)""" 

1900 BEEP 

1910 INPUT Ido 

1920 IF Ido>3 THEN Ido*3 

1930 IF ldo«0 THEN 3680 

19401 

19501 LOOP TO SET UPPER TUBE HEAT FLUX 

1960 IF Ido»l THEN 

1970 OUTPUT 709; "AK AF20 AL21 VR5" 

1980 PRINT "SET VOLT BOX TO MAIN" 

1990 BEEP 

2000 INPUT "ENTER DESIRED Qdp",Dqdp 

2010 PRINT USING "4X. ""DESIRED Qdp ACTUAL Qdp""" 
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2020 Err=1000 

2030 FOR 1=1 TO 2 

2040 OUTPUT 709; "AS SA" 

2050 Smn=0 
2060 FOR Ji=l TO 5 
2070 ENTER 709; E 
2080 Sum=Sura+E 
2090 NEXT Ji 

2100 IF 1=1 THEN Volt=Sum/5 
2110 IF 1=2 THEN Afflp=Sum/5 
2120 NEXT I 

2130 AjDp=ABS(Amp>» 1.9182) 

2140 Volt=ABS(Volt*25) 

2150 Aqdp=Volt^Amp/(PI*D2*L) 

2160 IF ABS(Aqdp-Dqdp)>Err THEN 

2170 IF Aqdp>Dqdp THEN 

2180 BEEP 4000,. 2 

2190 BEEP 4000,. 2 

2200 BEEP 4000,. 2 

2210 ELSE 

2220 BEEP 250,. 2 

2230 BEEP 250, .2 

2240 BEEP 250,. 2 

2250 END IF 

2260 PRINT USING "4X,MZ. 3DE, 2X,M2. 3DE" ;Dqdp,Aqdp 

2270 WAIT 2 

2280 GOTO 2030 

2290 ELSE 

2300 BEEP 

2310 PRINT USING "4X,MZ . 3DE , 2X,MZ . 3DE" ; Dqdp, Aqdp 

2320 Err=500 

2330 WAIT 2 

2340 GOTO 2030 

2350 END IF 

2360 END IF 

23701 

23801 LOOP TO SET Tsat 
2390 IF Ido=2 THEN 
2400 IF Ikdt=l THEN 2440 
2410 Dtld=2.22 

24201 PRINT USING "4X,"" DTsat ATsat Rate Tv 

2430 I)cdt=l 

2440 Oldl=0 

2450 Old2=0 

2460 Nn=l 

2470 Nrs=Nn MOD 15 

2480 Nn=Nn-a 

2490 IF Nrs=l THEN 

2500 PRINT USING "4X,"" Taat Tldl Tld2 Tv 

2510 END IF 

2520 OUTPUT 709; "AR AF16 AL19 VR5" 

2530 FOR 1=1 TO 6 
2540 IF I>4 THEN 2700 
2550 Sujn=0 

2560 OUTPUT 709; "AS SA" 

2570 FOR Ji*l TO 20 
2580 ENTER 709;Eliq 
2590 Sum=Sujn+Eliq 
2600 NEXT Ji 
2610 Ellq=Sum/20 
2620 Tld=FNTvsv ( Bliq) 



2630 


IF 


1 = 1 


THEN 


Tldl=Tld 


2640 


IF 


1=2 


THEN 


Tld2=Tld 


2650 


IF 


1=3 


THEN 


Tv-Tld 


2660 


IF 


1 = 4 


THEN 


Tsump=Tld 


2670 


IF 


1=5 


THEN 


Tinlet»Tld 


2680 


IF 


1=6 


THEN 


Tout»Tld 



Rate""" 



Tsump""" 
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2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

31801 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

laso 



NEXT I 

Atld*(Tldl+Tld2)*.5 

IF ABS(Atld-Dtld)>.2 THEN 

IF Atld>Dtld THEN 

BEEP 4000,. 2 

BEEP 4000,. 2 

BEEP 4000,. 2 

ELSE 

BEEP 250, .2 
BEEP 250, .2 
BEEP 250, .2 
END IF 

Errl=Atld-01dl 

01dl=Atld 

Err2=Tv-Old2 

01d2=Tv 

IF Tldl>100. THEN 2880 

PRINT USING "4X,5(MDD.DD,2X)";Dtld,Tldl,Tld2,Tv,T8ump 
WAIT 2 
GOTO 2470 
ELSE 

IF ABS(Atld-Dtld)>.l THEN 

IF Atld>Dtld THEN 

BEEP 3000,. 2 

BEEP 3000,. 2 

ELSE 

BEEP 800,. 2 
BEEP 800,. 2 
END IF 

Errl=Atld-01dl 

01dl»Atld 

Err2=Tv-01d2 

01d2=Tv 

PRINT USING "4X,5{MDDD.DD,2X)";Dtld,Tldl,Tld2,Tv,T3ump 

WAIT 2 

GOTO 2470 

ELSE 

BEEP 

Errl=Atld-01dl 

01dl=Atld 

Err2=Tv-01d2 

01d2=Tv 

PRINT USING "4X,5(MDDD.DD,2X)";Dtld,Tldl,Tld2,Tv,T3ump 

WAIT 2 

GOTO 2470 

END IF 

END IF 

END IF 

i 

LOOP TO SET LOWER TUBE HEAT FLUX 
IF Ido=3 THEN 

PRINT " SET VOLT BOX TO AUX" 

OUTPUT 709; "AR AF20 AL22 VR5" 

BEEP 

INPUT "ENTER DESIRED AuQdp" , Duxqdp 

PRINT USING "2X, ""DESIRED AuxQdp ACTUAL AuxQdp""" 

Err»1000 

FOR 1=1 TO 3 

OUTPUT 709; "AS 3A" 

Suio=0 

FOR Jl»l TO 5 
ENTER 709; E 
Sua=SiimtS 
NEXT J1 

IF I»1 THEN Volt=Sum/5 
IF 1=3 THEN Amp=Sum/5 
NE^ I 
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3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

36101 

3620 

3630 

3640 

3650 

36601 

3670! 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 



Ainp=AflS ( Amp* 1 . 9 182 > 

Volt=ABS(Volt*25) 

Auxqdp=Volt *Amp/ ( PI *D2 *L) 

IF ABS( Auxqdp-Duxqdp) >Err THEN 

IF Auxqdp>Duxqdp THEN 

BEEP 4000,-2 

BEEP 4000,-2 

BEEP 4000,-2 

ELSE 

BEEP 250,-2 
BEEP 250,-2 
BEEP 250, -2 
END IF 

PRINT USING "4X,MZ-3DE,2X,MZ-3DE";Duxqdp,Auxqdp 

WAIT 2 

GOTO 3260 

ELSE 

BEEP 

PRINT USING "4X,MZ-3DE,2X,MZ-3DE";Duxqdp,Auxqdp 

Err=500 

WAIT 2 

GOTO 3260 

END IF 

GOTO 1850 

END IF 

ERROR TRAP FOR Ido OUT OF BOUNDS 

IF Ido>3 THEN 

BEEP 

GOTO 1850 
END IF 



TAKE DATA LOOP 
IF Ikol=l THEN 3730 
BEEP 
Ikol=l 

INPUT ’’ENTER BULK OIL %”,Bop 
Ikol=l 

OUTPUT 709; "AR AFOO AL19 VR5” 
Ntc*20 



FOR 1=1 TO Ntc 
OUTPUT 709; "AS SA" 

Sum*0 

FOR Ji=l TO 20 
ENTER 709;E 
Sum=Sum+E 

IF I=(17-Nsub) OR I=(18-Nsub) THEN Et(Jl-l)=B 

NEXT J1 

Kdl=0 

IF I=(17-Nsub) OR I=(18-Nsub) THEN 

EavQ=Sum/20 

Sunj=0- 

FOR Jk=0 TO 19 

IF ABS(Et(Jk)-Eavo)<5-0B-6 THEN 
Suin=Sura+Et( Jk) 

ELSE 

Kdl=Kdl+l 
END IF 
NEXT Jk 

IF I=(17-Nsub) OR I=(18-Nsub) THEN PRINT USING 

IF Kdl>20 THEN 

BEEP 

BEEP 

PRINT USING "4X,""Too much scattering in data 
GOTO 1840 
END IP 
END IP 

Bmf(I)=Sum/{20-Kdl) 



4X,"’’Kdl = "",DD";Kdl 



repeat data set""" 
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4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

42901 

4300t 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

469Q 



NEXT I 
Coun^O. 

OUTPUT 709; '*AR AP20 AL22 VR5" 

FOR 1=1 TO 3 
OUTPUT 709; "AS SA" 

Suffl=0 

FOR Ji=l TO 5 
ENTER 709 ;E 
Sum=SumtE 
NEXT Ji 

IF Coun=0. THEN 
IF 1=1 THEN Vr=3um/5 
IF 1=2 THEN Ir=Sum/5 
ELSE 

IF 1=1 THEN Avr=Sum/5 
IF 1=3 THEN Air=Sum/5 
END IF 
NEXT I 

IF Coun=0 THEN 

PRINT "SHIFT VOLT BOX TO AUX" 

END IF 

IF Coun=l THEN GOTO 4290 

INPUT "TAKE AUX READINGS ( 1=YES) ?" ,Ccon 

IF Ccon>l THEN GOTO 4250 

Coun=Coun+Ccon 

GOTO 4060 

CONVERT emf'S TO TEMP, VOLT, CURRENT 

Twa=0 

Atwa=0 

FOR 1=1 TO Ntc 

IF IdtOO THEN 

IF I=Ldtcl OR I=Ldtc2 THEN 

T(I)=-99.99 

GOTO 4470 

END IF 

END IF 

IF AldtoO THEN 

IF I=Aldtcl OR I=Aldtc2 THEN 

T(I)=-99.99 

GOTO 4470 

END IF 

END IF 

T(I)=FNTvsv(Emf (I) ) 

NEXT I 

Twa=Twa/(4-Idtc) 

FOR 1=1 TO 8 

IF I=Ldtcl OR I=Ldtc2 THEN 

Twa=Twa 

ELSE 

Twa»Twa+T( I) 

END IF 
NEXT I 

Tw=Twa/(8-Idtc) 

FOR 1=9 TO 16 

IF I=AIdtcl OR I=Aldtc2 THEN 

Atwa=Atwa 

ELSE 

Atwa=Atwa+T( I ) 

END IF 
NEXT I 

Atw=Atwa/ ( 8-Aldtc ) 

Tld=T(17-Nsub) 

Tld2*T(18-Nsub) 

Tlda=(Tld+Tld2)’^.5 

Tv*T(19-Naub) 

Taump=T( 20-Nsub) 
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4700 

4710 

4720 

4730 

4740 

4750 

4760 

47701 

47001 

4790 

4000 

4010 

40201 

40301 

4040 

4050 

4060 

4070 

4000 

4090 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4900! 

49901 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5000 

5090 

5100 

51101 

5120! 

5130 

5140 

5150 

5160 

5170 

5100 

5190 

5200 

5210 

52201 

52301 

52401 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

53601 

53701 

53001 

5390 

5400 

5410 



Amp=ABS(Ir*1.9102) 

Volt*ABS(Vr) ^25 
Q=Volt*Ajnp 

Auanip=ABS( Air* 1 . 9102 ) 

Auvolt^ABS(Avr) <*25 
Auq=Auvolt*Auamp 
Kcu=Kcua( Itt) 

FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
Tw=Tw-Q^LOG(D2/Dl)/ (2*PI*Kcu*L) 

Atw=Atw-Auq’^LOG(D2/Dl) /(2’*PI*Kcu*L) 

Tsat=(Tlda+Tv)*.5 

COMPUTE THETHAB 
Th©tab=Tw-Tsat 
Athetab=Atw-Tsat 
IF Theta b<0 THEN 
BEEP 

INPUT ”UPPER TWALL<TSAT (0=CONTINUE, l=END)",Iev 
IF lev=0 THEN GOTO 3150 
IF I0V=1 THEN 1050 
END IF 

IF Athetab<0 THEN 
BEEP 

INPUT "LOWER TWALL<TSAT (0=CONTINUE, l=END)",Aiev 
IF Aiev=0 THEN GOTO 3150 
IF A1 qv=1 THEN 1850 
END IF 

COMPUTE VARIOUS PROPERTIES 

Tfllm=(Tw+Tsat)*.5 

Rho=FNRho(Tfilm) 

Mu=FNMu(Tfilm) 

K=FNK(Tfilm) 

Cp=FNCp(Tfilm) 

Beta=FNBeta ( Tf i Im) 

Hfg=FNHfg(Tsat) 

Ni=Mu/Rho 

Alpha=K/(Rho*Cp) 

Pr=Ni/Alpha 
Psat=FNPsat(Tsat ) 

COMPUTE VARIOUS PROPERTIES FOR AUX TUBE 
Atfilm=( Atw+Tsat) ♦ . 5 
Arho=FNRho(Atfilm) 

Ajtiu=FNMu( Atfilm) 

Ak=FNK(Atfilm) 

Acp=FNCp(Atfilm) 

Abeta=FNBeta ( Atf 1 Ira ) 

Anl=Arau/Arho 
Aalpha=Ak/ ( Arho* Acp ) 

Apr=Anl/Aalpha 

COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNBNHANCED END(S) 

Hbar=190 

Fe=(Hbar"P/(Kcu*A) )^.5*Lu 
Tanh=FNTanh(Fe) 

Th 0 ta=Thetab'^Tanh/Fe 

Xx=(9.81*B0ta*Thetab»Do"‘3*Tanh/(Fe*Ni^Alpha) )^. 166667 
Yy=( .559/Pr)^( 9/16) )^( 8/27) 

Hbarc=K/Do»( . 6+ . 307*Xx/Yy)^2 

IF ABS( (Hbar-Hbarc)/Hbarc)>.001 THEN 

Hbar=(Hbar+Hbarc)" .5 

GOTO 5260 

END IF 

COMPUTE NATURAL-CONVECTIVE HEAT -TRANSFER COEFFICIENT 

FOR UNENHANCED END(S) FOR AUX TUBE 

Ahbar=190 

F0= ( Ahbar "P/ ( Kcu"A ) ) " , 5 *Lu 
A tanh=FNTanh ( F© ) 
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5420 Atheta*Athetab*Atanh/Fe 

54 30 Axx=*(9. 81»Abeta»Athetab*Do"'3»Atanh/ (Fe^Anl^Adlpha) . 166667 

5440 Ayy=( lt( . 559/Apr)^( 9/16) )^( 8/27 ) 

5450 Ahbarc=Ak/Do* ( . 6+ . 387*Axx/Ayy) ^2 

5460 IF ABS( ( Ahbar-Ahbarc) /Ahbarc) > . 001 THEN 

5470 Ahbar=»( Ahbar+Ahbarc) * . 5 

5480 GOTO 5400 

5490 END IF 

55001 

5510! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 

5520 Ql=(Hbar*P*Kcu-A)^.5*Thetab*Tanb 
5530 Qc^Q-2*Ql 
5540 As=PI*D2*L 
5550 ! 



5560 1 COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) OF AUX TUBE 
5570 Aql=( Ahbar*P*Kcu*A) . 5*Athetab*Atanh 
5580 Aqc=Auq-2*Aql 
5590! 

56001 COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT 

5610 Qdp=Qc/As 

5620 HtubQ=Qdp/Thetab 

5630 Csf=(Cp*Thetab/Hfg)/(Qdp/(Mu*Hfg) *( . 014/ ( 9 . Bl^Rbo)^ . 5 )'' ( 1/3 . )*Pr^l.7) 

56401 

56501 COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT FOR AUX TUBE 

5660 Auqdp=Aqc/As 

5670 Ahtube=Auqdp/Athetab 

5680 Acsf=(Acp*Athotab/Hfg ) / ( Auqdp/ ( Aniu*Hfg) * ( . 014 / ( 9 . 81*Arho)'' . 5 ) ^ ( 1/3. ) *Apr-^l 
.7) 

56901 

57001 COMPUTE CHURCHILL-CHU CORRELATION FOR UPPER TUBE 

5710 Cc=Hbar*Thetab 

57201 

57301 RECORD TIME OP DATA TAKING 
5740 OUTPUT 709; ’’TD" 

5750 ENTER 709;Told$ 

57601 

57701 OUTPUT DATA TO PRINTER 
5780 PRINTER IS 701 
5790 PRINT 

5800 PRINT USING "10X,""Data Set Number = "" ,DDD, 2X, ""Bullc Oil % = , DD . D" ; J , B 



op 

5810 PRINT " 
5820 PRINT 



5830 


PRINT 


USING 


"10X,""TC No: 


8 


rt It tt 






5840 


PRINT 


USING 


"10X,""Temp : 


7),T(8) 






5850 


PRINT 


USING 


"10X,""TC No: 


16"**" 






5860 


PRINT 


USING 


"10X,""Temp : 


4),T(15),T(16) 




5870 


PRINT 






5880 


PRINT 


USING 


"lOX,"" Twa 


jjjpii"" 








5890 


PRINT 


USING 


"10X,3(MDD.DD 



TIME : " , TIME$ ( TIMEDATE ) 

1 2 

3( lX.MDD.DD] 
9 10 



11 



12 



13 



14 



15 



ld,Tld2 ,Tv,Psat,Tsat,Tsuxnp 



5900 PRINT USING "lOX,’”’ Vr Ir Cc”"** 

5910 PRINT USING "lOX, 2 ( MDDD .DD , 3X ) , 1X,M2. 3DE" ; Vr , Ir , Cc 

5920 PRINT USING "lOX,"" Thetab Htube Qdp Athetab Ahtuba AuQd 



p fV tv fV 

5930 PRINT USING "10X,MDD . 3D, 1X,MZ. 3DE, IX,MZ. 3DE, 1X,MDD. 3D, 1X,MZ. IDE, 1X,MZ . 3DE" 
; Xhatab, Htube , Qdp, Athetab , Ahtube , Auqdp 



5940 PRINT 
5950 BEEP 

5960 INPUT "OK TO STORE THIS DATA SET ( 1-Y, 0=N )7" ,01c 
5970 IP OX=l THEN J=J+1 
5980 IF Ok»l AND lm=0 THEN 

5990 OUTPUT @Filel;Bop,Told$,EiB£(*),Vr,Ir,Avr,Alr 
6000 END IF 

6010 IF Iuf-1 THEN OUTPUT @Ufile;Vw,Uo 
6020 IF Ire-1 THEN OUTPUT §RefHe;FmB,Rew 
6030 BEEP 
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6040 INPUT "WILL THERE BE ANOTHER RUN ( , 0=N ) , Go_on 
6050 Nrun*J 

6060 IF Go_on»0 THEN 6080 
6070 IF Go_on<>0 THEN Repeat 
6080 BEEP 
6090 PRINT 

6100 PRINT USING ”10X,""NOTE: data rune were stored in file 

l,D2_file$ 

6110 ASSIGN (BFilel TO * 

6120 OUTPUT (aFile2;Nrun-l 
6130 ASSIGN (apilel TO Dl_file$ 

6140 ENTER @Filel;Date$ ,Ldtcl,Ldtc2 , Itt 
6150 ENTER §Fllel; Aldtcl , Aldtc2 
6160 OUTPUT (aFile2;Date$,Ldtcl,Ldtc2, Itt 
6170 OUTPUT ^File2;Aldtcl,Aldtc2 
6180 FOR 1=1 TO Nrun-1 

6190 ENTER @Filel; Bop,Told$ , Emf ( , Vr , I r 

6200 ENTER @Filel; Avr , Air 

6210 OUTPUT ^File2;Bop,Told$,Emf(^) ,Vr,Ir 

6220 OUTPUT @File2; Avr, Air 

6230 NEXT I 

6240 ASSIGN 0Filel TO ♦ 

6250 PRINTER IS 1 
6260 PURGE "DUMMY” 

6270 PRINT 

6280 PRINT "DUMMY FILE PURGED" 

6290 SUBEND 
63001 

63101 CURVE FITS OF PROPERTY FUNCTIONS 
63201 

6330 DEF FNKcu(T) 

63401 OFHC COPPER 250 TO 300 K 

6350 Tk=T>273.l5 IC TO K 

6360 K=434-.112"Tk 

6370 RETURN K 

6380 FNEND 

63901 

6400 DEF FNMu(T) 

6410 COM /Iprop/ I ft 
6420 IF Ift=l THEN 
6430 1223 TO 373 K 

64401 170 TO 360 K CURVE FIT OF VISCOUSITY 

6450 Tk=T+273.15 IC TO K 

6460 Mu=EXP(-4.4636+(1011.47/Tk) ) ^l.OE-3 

6470 END IF 

6480 IF Ift=l THEN 

64901 223 TO 373 K 

6500 Mu=3.4593-4.26E-2*T>3.9485E-4^T^2-4.193E-6"T"3+2.0709E-8"T^4 

6510 Mu*Mu"1.0E-4 

6520 END IF 

6530 RETURN Mu 

6540 FNEND 

65501 

6560 DEF FNCp(T) 

6570 COM /Iprop/ I ft 

6580 IF lft=0 THEN 

65901 180 TO 400 K CURVE FIT OF Cp 

6600 Tk-T^'273.15 1C TO K 

6610 Cp=. 40188+ 1. 65007E-3"Tkf 1.5 1494E-6*Tk^2 -6.67 853B-10*Tk^3 

6620 Cp=Cp"1000 

6630 END IF 

6640 IF Ift=l THEN 

66501 223 TO 373 K CURVE FIT Cp(R124) 

6660 Cp=l. 0542+2. 1405E-3*T+1.0709E-5*T^2-6.4721B-8"T^3-1.4324E-9"T^4+ 4. 136E-11* 
T^5 

6670 Cp=Cp"1000 

6680 END IF 
6690 RETURN Cp 
6700 FNEND 
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67101 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

68101 

6820 

6830 

6840 

6850 

68601 

6870 

6880 

6890 

69001 

6910 

6920 

69301 

6940 

69501 

6960 

6970 

6980 

6990 

70001 

7010 

7020 

7030 

7040 

7050 

7060 

70701 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

72101 

7220 

7230 

7240 

7250 

7260 

7270 

72801 

7290 

7300 

7310 

7320 

7330 

7340 

73501 

7360 

7370 



DEF FNKho(T) 

COM /I prop/ I ft 
IF lft=0 THEN 
Tk=T+273.15 IC TO K 

X=l-(1.8^Tk/753.95) IK TO R 

Ro=36. 32+61. 146414*X^( 1/3) +16. 418015*X+ 17. 476838*X^. 5+1. 119828^X^2 
Ro=Ro/. 062428 
END IF 

IF Ift=l THEN 
223 TO 373 K (R124) 

Ro=1434.8-2.8619*T-6.7267E-3*T^2-7.2852E-5*T^3 
END IF 
RETURN Ro 
FNEND 

DEF FNK(T) 

COM /I prop/ Ift 

IF lft=0 THEN 

T<360 K WITH T IN C 

K=.071-.000261*T 

END IF 

RETURN K 

IF Ift=l THEN 

223 TO 373 K (R-124) 

K=7.5191E-2-3.5436E-4*T-1.9545E-7*T^2+3.1835E-9*T^3 
END IF 
RETURN K 
FNEND 

DEF FNTanh(X) 

P=EXP(X) 

Q=l/P 

Tanh=(P-Q)/(P+Q) 

RETURN Tanh 
FNEND 

DEF FNTvsv(V) 

COM /Cc/ C(7),Ical 
T=C(0) 

FOR 1=1 TO 7 
T=T+C(I)*V-I 
NEXT I 

IF Ical=l THEN 

T=T-6.7422934K-2+T^(9.0277043E-3-T*(-9.3259917E-5) ) 

ELSE 

T=T+8.626897E-2+T^(3.76199K-3-T*5.0689259E-5) 

END IF 
RETURN T 
FNEND 

DEF FNBeta(T) 

Rop=FNRho(T+ . 1) 

Rom=FNRho(T-.l) 

Beta=-2/ (Rop+'Rom) *(Rop-Rom) / . 2 

RETURN Beta 

FNEND 

DEF FNHfg(T) 

COM /I prop/ Ift 
IF lft=0 THEN 

Hfg»1.3741344K+5-T*( 3. 3094361E+2+T* 1.2165 143) 

END IF 

IF Ift=l THEN 

223 TO 373 K CURVE FIT (R124) 

Hfg=15 9.93-.4609*T-1.458E-3’*T^2-1.3715E-5*T^3 
Hfg=Hfg*1000 
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7380 2ND IF 
7390 RETURN Hfg 
7400 FNEND 
7410! 

7420 DEF FNPsat(Tc) 

7430 COM /Iprop/ Ift 
7440 IF lft=0 THEN 

74501 0 TO 80 deg F CURVE FIT OF ?sat 
7460 Tf='1.3*Tc*^32 

7470 ?a=5.945525^Tf^( . 15352082>Tf* ( 1 . 4840963E-3^Tf *9 . 615067 lB-6 ) ) 

7480 ?g*Pa-14.7 

7490 END IF 

7500 IF Ift^l THEN 

75101 223 TO 373 K CURVE FIT (R-124) 

7520 Pa-( 162.68^-5.946*Tc-^9.2081E-2*Tc"'2 + 6.9359E-4^Tc"‘3)*. 14504 
7530 Pg*Pa-14.7 

7540 END IF 

7550 IF Pg>0 THEN 1 +»PSIG,-=tn Hg 

7560 Psar=«Pg 
7570 ELSE 

7580 Psai:=Pg»29.92/14,7 

7590 END IF 
7600 RETURN Psat 
7610 FNEND 
7620 : 

7630 i 

7640 SUB Purg 
7650 3EEP 

7660 INPUT "ENTER FILS NAME TO 3E DELETED" ^FileS 
7670 PURGE File$ 

7680 GOTO 7650 
7690 SUBEND 
7700 : 

7710 : 

7720 SUB Select 
7730 COM /Idp/ Idp 
7740 3EEP 
7750 PRINTER IS 1 

7760 PRINT USING "2X, ""Select ootion;""" 

7770 PRINT USING "6X,"" 0 TAKE DATA""" 

7780 PRINT USING "SX,"" 1 PURGE FILES""" 

7790 INPUT Idp 
7800 IF Idp-0 THEN GALL Main 
7810 IP Idp«l THEN CALL Purg 
7820 SUBEND 
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APPENDIX C. REPRESENTATIVE DATA SET 



Jdta •* 7 SeD 1993 

NOTE: Program nane : ORPJY 

New rile na/ne: O0907SC2 

UPPER *U8E TC are defective at LOCATIONS; 7 3 



t 


OUER TUBE TC : 


5 defective at L 


.OCATION; 9 




1 


uQea: Uieland 


Hard 3 incn 






Oat a Set Nunoer * 1 


3uik Oil •: » 


0.0 




TIME: 


13: 15:35 








TO No: 


1 2 


3 4 


5 o 


7 3 


T enp : 


3.35 3.73 


3.75 3. So 


3.37 3.39 - 


99.99 -99.99 


TC No: 


9 10 


11 12 


13 14 


15 IS 


T eng : 


-99.99 20.20 


13. 9S 20. 4S 1 


9.54 20.49 


13.53 17.30 


T wa 


ATwa Tiiqd 


TiiqdZ Tvapr 


Psat Taat 


T 3unp 


3.79 


19. 3S 2.20 


2.22 2.11 


10.32 2.15 


-10,3 


Ur 


Ir 


Cc 






.34 


. 17 


2.333£>02 






Thetao 


Htube 


Qdp Athetab Ahtube 


AuQdp 


1 .S2S 


3.403E+02 5 


. 5425+02 17.194 


1 .3325-02 


3. 1505-03 



Oata Set 


Nup^ber * 2 


3uik Oil r. 0.0 








TIME: 


13: 


14:33 










TC No: 


1 


2 


3 4 5 


5 


7 


9 


Temp : 


4.37 


4.91 


4.32 4.32 4.9S 


5.10 


-99.99 - 


99.39 


TC No: 


9 


10 


n 12 13 


1 4 


15 


16 


"^enp : - 


99.99 


20.49 


19.22 20.57 20.07 


20.32 


13.33 


17.53 


Tua 


ATwa 


T1 iqd 


TliQd2 Tvapr Psat 


Tsat 


T sung 




4.94 


19.55 


2.^9 


2.36 2.13 10.90 


2 . 2 


5 -10.3 




Ur 




Ir 


Cc 








1.13 




.23 


4 . 530E+02 








Thetao 


Htube 


Qdp Athetab Ahtube 


AuQdp 




2.533 


3.3G3E+02 1. 


.040E+03 17.397 1.793E-02 


3.120E+03 



Oata Set 


Number =• 3 


Bulk Oii X * 


0.0 




TIME: 


13:16:14 








TC No: 


1 2 


3 4 


5 6 


7 9 


Tenp : 


7.38 7.12 


7.12 5.32 


7.19 7.43 - 


99.99 -99.99 


TC No: 


9 10 


11 12 


13 14 


15 16 


Tenp : - 


99.99 20.08 


13.76 20.32 


19.41 20.16 


13.35 16.96 


Tua 


ATwa Tliqd 


Tliqd2 Tvapr 


Psat Tsat 


T sunp 


7.17 


19.12 2.07 


2.13 2.02 


10.73 2.06 


-10.3 


Ur 


Ir 


Cc 






1 .49 


.31 


1 . 02SE+03 






Thetao 


Htube 


Qop ftthetab rthtuoa 


AuQdp 


5.104 


3.4105^02 1. 


,740E-.03 17.357 I.380E+02 


3 . I 7ZE+03 
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Qata Set 


NunOer * 4. 


3ulk Oil I - 0.0 






TIi'lE: 


13: 


17;a3 








TC No: 


I 


z 


3 t S 


5 


7 3 


T enp : 


1 3. 4.0 


10.10 


10.12 9.54 10.01 


1 

m 


qq^qg _gg qg 


TC No: 


9 


10 


11 12 13 


1 4 


15 15 


Tanp : - 


99.99 


20.02 


13.77 20.12 19.51 


20.29 


13.32 17.29 


Twa 




T 1 iqd 


TLiqd2 Tvapr Psat Taat 


T 3unp 


13.11 


19.13 


1 .95 


2.15 1.92 10. 


56 1.99 


-10.7 


^Jr 




Ir 


Cc 






1 .37 




.39 


1 , 855£'*-03 






Thetao 


Ht ube 


Qdo Athetao AhtuPe 


AuOdp 


3.121 


3.3196^02 2. 


,595£'"03 17.139 1. 


3346-02 


3. 1525-03 



Data Set 


Nunber * 5 


3uU Oil r. » 0.0 




TL^IE: 


13:13:54 








TC No: 


1 2 


3 


4 5 5 


7 3 


Tenp : 


14.37 14.17 


14.12 


13.47 13.75 14.60 - 


99.99 -99.39 


TC No: 


9 10 


1 1 


12 13 14 


15 16 


Tenp ; - 


99.99 20.49 


19.23 


20.42 20.02 20.52 


13.47 17.57 


T’wja 


ATwa Tliqd 


Tl iqd2 


Tvapr Psat Taat 


T 3unp 


14.05 


19.52 2.01 


2.25 


2.04 10.76 2.09 


-10.3 


Vr 


Ir 


Cc 






2.36 


.50 


3.0465- 


03 




That ao 


Ht ube 


Qdo 


^thetab Hhtube 


AuOdp 


11.360 


3.530£t02 4, 


.3425^03 


17.429 1.7306+02 


3. 1025-03 



Data Set 


NunOer 


' * 3 


9ulk Oil ^ - 0 


.0 




TIfl6: 


13:20 


i:03 








TC No: 


1 


2 


3 4 


5 3 


7 3 


Tero : 


17.39 


17.40 


17.30 16.51 16 


.34 17.33 " 


99.99 -^9.99 


TC No: 


9 


10 


1 I 12 


13 14 


15 IS 


Tenp : - 


cn 

m 

cn 

cn 


20.21 


13.38 20.32 19 


.39 20.48 


18.47 17.08 


Tua 


i^Twa 


Tl iqd 


Tliqd2 Tvapr 


Peat Taat 


T 3uno 


17.16 


19.27 


2.04 


2.32 2.09 


10.30 2.13 


-10.7 


\Jr 




Ir 


Cc 






2,53 


. 


55 


4 . 0846+03 






Thetao 


Htuba 


Qdp ^thetab 


f^htube 


AuOdp 


15.024 


3.3796+02 5 


.0776+03 17.139 


1 .3346+02 


3.1436+03 
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7 



3uU 3il 



0.3 



Qdta Set Numoer = 
TIME: 13:21:15 



1 

o 

o 


1 


2 


3 >i 5 3 


7 3 


emo : 


19. -i3 


19.13 


13,95 13.35 13.12 19.37 - 


99.99 “99.99 


TC No: 


9 


10 


11 12 13 U 


15 IS 


Tema : 


“99.99 


19.99 


18.74 19.94 19.50 20.11 


13.38 17.37 


TuJd 


P»TuJd 


rUqd 


riiqaZ 7vaqr Psdt Tsat 


Taumq 


18,79 


19.33 


1 .39 


2.21 I. 96 10.58 2.30 


“10.7 


Ur 




Ir 


Cc 




Z.30 


, 


30 


4.7075+03 




Thetao 


Htube 


Qdp Athetab Ahtube 


AuQdo 


16.783 


3. SIZE 


-02 6 


.3615-^03 17.327 1.348E+02 


3 . 1 46E-03 



Data Set 


Number » 8 


Bulk Oil r. - 0.0 








rr.NE: 


I3;2Z: 1 1 












TC No: 


1 2 


3 


4 5 6 


7 




3 


Temp : 


22.76 22.21 


21 .96 


20.36 21 .04 22.51 - 


99.99 - 


99 


.99 


TC No: 


9 10 


1 1 


12 13 14 


15 




IS 


T erp : - 


99.99 19.75 


18.48 


19.79 19.35 20.35 


18.04 


16 


i/i 

CD 


Tua 


/^Tua Tliqd 


TI iqd2 


Tvapr Paat Taat 


T sump 






21 .34 


18.87 1.31 


2.19 


1.90 10.63 1.95 


-10.7 






Ur 


Ir 


Cc 










3.37 


.65 


5 . 3555+ 


03 








Thetab 


Htube 


Qdp 


Athetab Ahtube 


AuCdp 






19.389 


3.552E^02 7. 


. 263E+03 


16.325 I.364E-02 


3.1555* 


03 





Data Sat 


Number » 9 


Bulk Oil X - 0.0 




TINE: 


13: 


23:55 






TC No: 


1 


2 


3 4 S 6 


7 8 


Temp : 


2S.31 


26.51 


26.32 25.12 24.95 26.71 


-99.99 -^9.99 


TC No.- 


9 


10 


11 12 13 14 


15 IS 


Temp : - 


99.99 


20.53 


19.25 20.45 20.10 20.75 


18.59 17.71 


Tya 


A Tuid 


Tliqd 


TliqdZ Tvapr Psat Taat 


T sump 


26.03 


19.30 


2.14 


2.47 2.23 10.32 2.2 


7 -10.7 


Ur 




Ir 


Cc 




3.29 




.70 


7.3S9E+03 




Thatab 


Htube 


Qdp Athatab Ahtuba 


AuQdp 


23.757 


5.331E+02 3. 


,05SE+03 17.523 1.785E+02 


3.094E+03 
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Data Set Nuinoer ^ !3 9uik Oil *1 » 0.0 



TIME: 


13:27 


’:26 










rC No: 


I 




3 


4 


5 5 


7 3 


Teno : 


10.92 


1 1 .50 


1 1 .27 


11.74 11 


.29 10.35 - 


99.99 -99.99 


TC No: 


9 


10 


1 1 


12 


13 14 


15 16 


Teno : ' 


-99.99 


17.35 


15.33 


17.30 17 


.07 13.05 


15.48 14.36 


Tua 


Aiwa 


T 1 iqd 


71 iqd2 


Tvapr 


Psat Tsat 


7 3urp 


11.16 


16.34 


2. 23 


2.25 


2.09 


10.33 2.17 


-10.5 


Ur 




Ir 


Cc 








3.79 




31 


2. 1 I4£t03 






Thetab 


Ht ube 


OdD 


Athotao 


Ahtube 


AuOdp 


3.991 


1 .510E 


^03 I 


.357E+04 


14.471 


2.4S5E^02 


3.553Er03 



Data Sot 


Nunber » 1 1 


3ulk Oil 0.3 






7I.'1£: 


13:29:22 








o 

o 

1 — 


1 2 


3 4 5 


6 


7 3 


7erp : 


11.77 12.52 


12.23 12.73 12.40 


1 1 .91 


-99.99 -99.99 


7C No: 


9 10 


11 12 13 


1 4 


IS 16 


7enp : - 


99.99 16.90 


15.33 16.93 16.50 


17.31 


16.36 13.76 


7ya 


A7ua 7liqd 


7liqd2 7vapr Psat 


7aat 


7 3urp 


12.14 


16.21 2.23 


2.22 2.01 10.73 2.1 


2 -10.5 


Ur 


Ir 


Cc 






4.43 


.96 


2.430E-03 






7hetab 


Htube 


Qdp Athetab Ahtube 


AuOdp 


10.027 


I.924E-03 1 


.9295^04 14.092 4.2 


13E^02 


5 . 944£*<*03 



Data Set 


Nurber 


- 12 


3ulk Oil 7. 


0.0 






7IflE: 


13:30 


:20 










7C No: 


1 


2 


3 4 


5 


6 


7 - 3 


7erp : 


13.30 


14.46 


14.14 14.73 


14.52 


13.94 - 


■99.99 -99.99 


7C No: 


9 


10 


1 1 12 


13 


14 


15 16 


7efip ; - 


'99.99 


16.21 


14.70 16.26 


15.95 


17.07 


15.51 13.21 


7ua 


ATua 


71 iqd 


71iqdZ 7vapr 


Psat 


7sat 


7 


1 4.07 


15.53 


Z. I 1 


2.14 1.93 


10.70 


2.03 


-10.5 


Ur 




Ir 


Cc 








5.33 


1 . 


25 


3.072E403 








7hotab 


Htube 


Qdp Athetao Ahtube 


AuOdp 


12.039 


2.792E+03 3. 


.jSIE^-04 13.502 2.7S5E+02 


3.720E+03 
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Data Sat 
TIME: 


Numoer =» 1 3 

13:31 :4j 


3u;k Oil -C - 0.0 




rC Mo: 


I : 


3 4 9 5 


7 3 


^ arp : 


16.12 15. ■’5 


15.45 17.02 15.97 16.17 


-99.99 -99.99 


rC No: 


9 10 


11 12 13 14 


15 15 


^emp : - 


99.39 15. “’S 


M.44 15.97 15.71 15.93 


15.49 13.23 


^wa 


(4Tuja TliQd 


71iqd2 Tvapr Psat Taat 


T 9una 


!S.:4 


15.34 2.21 


2.32 2.16 10.37 2.2 


1 -10.5 




Ir 


Co 




7.71 


1 .54 


3.739Ei'03 




Thetad 


Ht uae 


Qap f^thetad rthtaba 


^uOdo 


U. J30 


4.162E-03 5. 


.333E-04 13.125 2.369E+02 


3.756E+03 



Data Set 


NuDOer » 1 4 


9ulk Oil X - 0.0 






TIME: 


13: 


32:43 








rC No: 


1 


-7 


3 4 5 


5 


7 9 


7 emo : 


17.77 


18.27 


17.93 13.48 18.53 


17.70 - 


99.99 -99.99 


rC No: 


9 


10 


11 12 13 


1 4 


15 16 


Tsno : - 


99.99 


15.53 


14.22 15.55 15.39 


15.51 


14.99 12.90 


T ya 


f^Twa 


T I iqd 


TIiqd2 Tvapr Psat 


Tsat 


T sump 


17.34 


15.03 


2.29 


2.36 2.21 10.92 


2.27 


-10.5 


Vr 




Ir 


Cc 






9.23 




1 .96 


4.207E+03 






Thetaa 


Ht ube 


Qdo i^thetab Phtube 


AuOdp 


15.375 


5.453E+03 9, 


,33i£-04 12.763 2.37 


2E-02 


3.7925^03 



Data Set 


Number 


- 15 


3uik Oil X - 0.0 






riflE: 


13:34 


:20 








TC No.- 


1 


7 


3 4 5 


5 


7 - 3 


Temp : 


16.36 


17.52 


17.20 17.74 17.59 


16.90 - 


99.99 -99.99 


TC No: 


9 


10 


11 12 13 


1 4 


15 15 


T emp : - 


99.39 


15.59 


14.13 15.58 15.39 


16.54 


15.18 12.94 


Tya 


PTya 


TI iqd 


TliqdZ Tvapr Psat 


Tsat 


T sumo 


15.95 


15.05 


2.40 


2.54 2.31 11.03 


2.39 


-10.4 


Vr 




Ir 


Cc 






3.23 


1 . 


75 


3.925E+03 






Thetad 


Htubo 


Qdp Pthetab Phtube 


PuQdo 


14.563 


4-.572E+33 S. 


.5536^04 12.562 2.987E+02 


3.782EH-03 
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Data Sat 


Nunoer = IS 


9uU Oil X - 0.0 


TIME: 


13: 


35:45 




-1 

o 

o 


1 


7 


345379 


T ewQ ; 


15.51 


16. :9 


15.38 15.59 16.32 15.30 -99.39 -39.99 


TC Mo: 


q 


10 


II 12 13 14 15 16 


Tano : - 


■99.99 


15.49 


14.01 15.51 15.31 16.33 U.36 12.72 


T uja 


r^Tuia 


Tl iqd 


TLiqdZ Tvaor ?sat Tsat Tauno 


15.30 


U.36 


2.33 


2.35 2.25 10.34 2.29 -!0.4 






Ir 


Cc 


5.38 




1 .47 


3 . 5o0E't*03 


Thatab 


Htube 


Qdo i^thetab Ahtube AuQdp 


13.503 


3.43IE-M33 1 


.534E+’04 12.553 3.038E'^02 3.318E*03 



Data Set 


Number 


- 17 


Bulk Oil 5 : - 0.0 




TIME: 


13:36 


:S1 






TC No: 


1 




3 4 5 5 


7 3 


Temo : 


13.76 


14.46 


14. 18 14.31 14.21 13.35 - 


99.39 -99.39 


TC No: 


9 


10 


11 12 13 14 


15 16 


T eno : - 


■99.99 


15.21 


13.75 15.17 15.10 16.15 


14.50 12.42 


Tua 


ATua 


T i iqd 


Tliqd2 Tvaqr Paat Taat 


T aumq 


14.02 


14.50 


2.25 


2.29 2.08 10.33 2.17 


-10.2 


\)r 




Ir 


Cc 




5.28 


1 . 


13 


3.010E+03 




Thetab 


Htube 


QdD ^thetab r^htube 


AuQdo 


11.348 


2.271E^03 2. 


, 590E”‘*04 12.427 3.3955*02 


3.347E-03 



Data Set 


Number * 18 


Bulk Oil X - 0.3 








TIME: 


13: 


38:09 












TC No: 


1 


2 


3 


4 5 


6 


7 


8 


Tenq : 


11.36 


12.54 


12.42 


12.38 12.20 


11.69 - 


39.39 - 


99.99 


TC No: 


9 


10 


1 1 


12 13 


14 


15 


16 


Temp : - 


39.99 


15.29 


13.32 


15.29 15.25 


16.45 


15.37 


12.51 


Tua 


ATua 


Tl iqd 


TUqd2 


Tvapr Paat 


Tsat 


Taumo 




12.21 


14.30 


2.28 


2.24 


2.06 10.32 


2.16 


-10. I 








Ir 


Cc 










4.36 




.87 


2.437E-^03 








Thetab 


Htube 


Qdp 


ftthetatj Ahtube 


^uOdp 




10.048 


1.5S6£-^03 1 


.5635+04 


12.333 3.305c-^02 


3. 7975^03 
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Odta Sat 


Nu no ar ® I 9 


9ulk Oii X = 0.0 






TL'IE: 


13:59:23 








TC No: 


1 2 


3 4 5 


5 7 


9 


"i* anp : 


9.20 9.53 


9.45 9.30 9.35 


9.39 -99.99 - 


39.99 


TC No: 


9 10 


n 12 13 


14 15 


1 3 


Tsno : - 


99.99 15.53 


14.21 15.4S 15.53 1 


15.38 15.32 


13.15 


T ya 


(9Tya Tilqd 


TiiqoZ Tvapr Psat 


Taat Tsunp 




9.50 


15.13 2.28 


2.35 2.17 10.90 


2.24 -10.1 




Ur 


Ir 


Cc 






2.72 


.53 


1 .5095^03 






Thetab 


Htuba 


Qdo Athatao Ahtuoa ^^uQdp 




7.251 


9.2455-02 5. 


,7125+03 12.390 2.3335-02 3.596c- 


03 



Data Sat 


Nunbar =* 20 


3uU Oil r. - 0.0 






TIME: 


13:40: 18 










TC No: 


1 2 


3 


4 5 


5 


7 3 


Tanp : 


7.12 7.13 


7.03 


7.10 7.79 


7.23 - 


99.99 -99.99 


TC No: 


9 10 


1 1 


12 13 


1 4 


IS 16 


Tanp : - 


99.99 15.55 


14.29 


15.45 15.33 


16.57 


15.31 13.21 


T ua 


ATwa Tliqd 


Tl iqd2 


Tvapr Psat 


Taat 


T 3unp 


7.21 


15.15 2.25 


2.32 


2.22 10.91 


2.25 


-10.1 


Ur 


Ir 


Cc 








2.04 


• .43 


9.3736+ 


02 






Thetab 


Htube 


Qdp 


Athatab i^htuoe 


(^uOdD 


4.955 


7. 4996-02 3, 


.7176+03 


12.900 2.9226-02 


3.7536-05 



Data Set 


Nunber 


• • 21 


Bulk Oil *: - 


0.0 






TIME: 


13:41 


:00 












TC No: 


1 


2 


3 


4 


5 


6 


7 8 


T enp : 


5.17 


6.12 


6.04 


6.02 


5.77 


6.33 - 


■99.99 -99.99 


TC No: 


9 


10 


1 I 


12 


13 


14 


15 16 


Temp ; - 


99.99 


15.78 


14.52 


15.53 1 


5.39 


16.76 


15.50 13.48 


Tua 


i^T ya 


T 1 iqd 


Tliqd2 


Tvapr 


Paat 


Taat 


T aunp 


3.23 


15.34 


2.29 


2.36 


2.25 


10.94 


2.29 


-10.0 


Ur 




Ir 


Cc 










1 .74 




37 


7.3676+02 








Thetab 


Htube 


Qdp 


Athatab /^htuba 


AuQdp 


3.33S 


6.8016+02 2 


.5776+03 


13.053 


2.3686+02 


3 . 7446+03 
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Data 3et 


Nurber 


a T 


3uik Oil 


1 » 0.0 






7inE; 


13:42 


: 13 










^o: 


1 


: 


3 


4 5 


9 


7 a 


T emp : 


4.33 


4 . 40 


4,39 4 


.37 4.33 


4.51 - 


99.99 -99.99 


o 

o 


9 


10 


1 1 


12 13 


1 4 


15 15 


Teno : “ 


■99.99 


15.38 


14.70 15 


.93 19.37 


19.95 


15.53 13.55 


T uia 


ATwa 


TL ’ qd 


TLiqd2 


Tvapr Psat 


Tsat 


7 3urp 


-t.i? 


15.51 


2.45 


2.51 


2.35 11.09 


2.41 


-10.1 


Ur 




Ir 


Cc 








1 .33 




27 


3.2255^02 








Thetao 


Htube 


Qaq 


Athetao f^htube 


.-^uOdp 


2 . 353 


7.540E>32 1. 


.572E+03 


13.395 2.393£>02 


3 . 7555+03 



Data Set 


Nurber =* 23 


3ulk Oil r. - 0.0 








ri.lE: 


13: 


43:07 












TC No: 


1 


2 


3 


4 5 


9 


7 


3 


— 1 

3 

O 


3.53 


3.53 


3.52 


3.47 3.33 


3.30 - 


•99.39 - 


99.99 


TC No: 


9 


10 


1 1 


12 13 


14 


15 


19 


Terp : - 


99.39 


15.38 


U.54 


15.38 15.98 


15.39 


15.54 


13.42 


Twa 


ATwa 


TI iqd 


TIiqd2 


Tvapr Psat 


Tsat 


T sunp 




3.53 


15.42 


2.31 


2.37 


2.19 10.92 


2.2? 


-10.0 




Ur 




Ir 


Cc 










.99 




.20 


1 .3255- 


02 








Thetao 


Htube 


Qdp 


r^thetao i^htube 


^uQdP 




1 .317 


3.52' 


9E*^02 3. 


597E-02 


13.153 2.308E-02 


3.593E- 


03 



Data Set 


Number 


- 24 


3ulk Oil 


r. » 0.0 




TIME: 


13:44 


:32 








TC No: 


1 


2 


3 


4 5 5 


7 _ 3 


Terp : 


Ui 

cn 


2.94 


2.95 2 


.33 3.12 2.99 - 


99.39 -99.39 


TC No: 


9 


10 


1 1 


12 13 14 


15 15 


Terp : - 


99.39 


15.99 


14.53 15 


.37 15.08 17.01 


15.93 13.50 


Tua 


hTua 


T 1 iqd 


Tl iqd2 


Tvapr Psat Tsat 


T sump 


2.38 


15.52 


2.45 


2.48 


2.34 11.05 2.40 


-10.0 


Ur 




Ir 


Cc 






.70 




1 4 


5.344E+0I 






Thetao 


Htube 


Qdp 


Athetao AhtuOe 


^uOdo 


.574 


7.393E 


^02 4. 


.Z41E-K02 


13.120 2.347E+02 


3.7355+03 


NOTE: 24 


data runs '^ere stored 


in file O0907SC2 
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APPEa^DIX D. SAMPLE CALCULATIONS^ 



Data file D0921ASC4 was used for sample calculations. The 
saturation temperature was 2.2 °C with a heat flux of 59,130 
W/m^ on the smooth upper test tube. 

A. SMOOTH TEST TUBE DIMENSIONS 

Do = 0.015875 m 
Di = 0.0127 m 
D1 = 0.012446 m 
L = 0.2032 m 
Lu = 0.0762 m 

B . MEASURED PARAMETERS 

V = Vr-25 = (7.76 V) (25) = 194.0 V 
I = Ir-1.9182 = (1.66 A) (1.9182) = 3.18 A 
T1 = 19.44 °C T2 = 20.03 °C 

T3 = 19.75 °C T4 = 19.63 °C 

T5 = 19.46 °C T6 = 19.38 °C 

T7 and T8 not used, defective 



^This analysis is essentially identical to the analysis 
done by Perry (Ref. 2] and therefore his procedure is largely 
reproduced. 
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Tsat = 2.20 °C 
kc = 344 W/m-K 

C. OUTER WALL TEMPERATURE OF THE UPPER TEST TUBE 

1. Surface Area of Smooth Tube 

P = Ji'Do = n(0. 015875 m) = 0.04987 m 

Ac = (Do^-Di^) /4 = ji((0. 015875 m) ^- ( 0 . 0127 ) ^ ) /4 

Ac = 7.13 X 10-5 

2 . Heater Power 

Qh = V-I = (194.0 V) (3.18 A) = 616.92 W 

3 . Average Inside Tube Temperature 

Tavg = Sr Tn/m = Ei® Tn/6 = (T1+T2+T3+T4+T5+T6) /6 
Tavg = (19.44+20.03+19.75+19.63+19.46+19.38) /6 
= 19.62 °C 

4 . Outside Wall Tube Tenperature 

Two is calculated using Fourier's Conduction Law, 
assuming uniform radial conduction and using a known Tavg: 

Two = Tavg - (Qh ( In (Do/Di ))/ (2n -L • kc) } 

Two = 19. 62-{616.92(ln(0. 015875/0. 01245) )/ 

(2ji-0. 203-344 ) } 

Two = 19.28 °C 

5. Tube Wall Superheat 

0 = Two - Tsat = 19.28 - 2.20 = 17.08 °C 



156 



D. R-124 PROPERTIES AT FILM TEMPERATURE 



The thermophysical properties of R-124 were computed at 
the film temperature using the REFPROP equations provided by 
Bertsch [Ref. 1] : 

1 . Film Ten5>erature 

Tf = (Two + Tsat)/2 = (19.28 + 2.20)/2 = 10.74 °C 

2 . Dynamic Viscosity 

p = 3.4593-(0.0426)Tf+(3.9485xl0-')Tf2-(4.193xl0-^)TfH 
(2.0709xl0-®)Tf^ 

p = (3.4593- (0.0426) 10.74+ (3. 9485x10-'') (10. 74)2- 
(4.193x10-^) (10. 74)2+(2. 0709x10-®) (10.74)^x10-' 
p = 3.0424x10-' N-s/m2 

3. Density 

p = 1434.8-(2.8619)Tf-(6.7267xl0-®)Tf2-(7.2852xl0-®)Tf® 
p = 1434.8- (2.8619) 10.74- (6.7267x10-2) (10. 74)2- 
(7.2852x10-') (10.74)2 
p = 1403.2 kg/m® 

4. Kinematic Viscosity 

V = p/p = 3.0422x10-' /1403.2 = 2.17X10'® m/s 

5 . Thermal Conductivity 

k = 7.5191xl0-2-(3.5436xl0-')Tf-(1.9545xl0-')Tf2 + 
(3.1835xlO-®)Tf2 

k = 7. 5191x10-2-(3. 5436x10-') (10. 74)-(l. 9545x10-2) 
(10.74)2+ (3.1835x10-®) (10.74)® 
k = 7.137x10-2 W/m-K 
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6 . 



Specific Heat 

Cp = 1.0542+(2.1405xl0-^)Tf+(1.0709xl0'^)Tf'- 

(6.4721x10'®) -Tf®- (1.4324x10'®) TfH (4 . 136x10'" )Tf® 

Cp = (1.0542+ (2.1405x10'®) ( 10 . 74 ) + ( 1 . 0709x10'® ) (10.74)®- 
(6.4721x10'®) (10.74) ®- (1.4324x10'®) (10.74)' + 
(4.136x10'®®) (10.74)®}xl0® 

Cp = 1078.35 J/lcg-K 

7. Thermal Diffusivity 

a = lc/(p'Cp) = 7.137x10'®/ (1403.2 -1078.35) 

a = 4.717x10'® m®/s 

8. Voliimetric Thermal Expansion Coefficient 

P = -Ap/(p-ATf) = - (1402.89-1403.50) / (1403.2 -0.2) 

P = 2.174x10'® (1/K) 

9 . Prandtl Number 

Pr = v/a = 2. 17xl0'®/4. 717x10'® = 4.60 

E. HEAT FLUX THROUGH NON-BOILING TUBE LENGTHS 

Using the Churchill-Chu correlation [Ref. 9] for a smooth 
cylinder in natural convection, the average natural convection 
heat transfer coefficient of one non-boiling tube length is: 

h = k/Do(0.6+0.387[r/(l+(0.559/Pr)®"®)®^®®] }® 

where F= { [g • P -Do® -0 • tanh (m-Lu) ] / (v • a -Lu *m) }®^® 
m = { (h-P) / (Icc-Ac) }®^® 
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The program DRPJY computes the value of h by assuming a 
starting value of 190 W/m^*K and iterating to the final value. 
The resulting values of h and m are: 

h = 281.93 W/m^-K 

m = 23.94 (1/m) 

therefore, 

Qf = (h 'p ’kc ’Ac) *0 *tanh (m*Lu) = 9.52 W 

F. HEAT FLUX THROUGH BOILING TUBE LENGTH 

Q = Qh - 2’Qf = 616.92 - 2(9.52) = 597.88 W 

Ab = ji’Do’L = 11(0. 015875) (.2032) = 1.013x10-^ m^ 
q" = Q/Ab = 597.88/1.013x10-2 = 59, 021 W/m^ 
h = q"/0 = 59021/17.08 = 3456 W/m^-K 

The following results were calculated by the data 
acquisition and reduction program DRPJY: 

q" = 59, 130 W/m2 

h = 3463 W/m2’K 

0 = 17.08 °C 
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APPENDIX E. UNCERTAINTY ANALYSIS 



4 - 

The same data run used for the sample calculations 
(Appendix D) was used for the uncertainty analysis, thus all 
the measured and most of the calculated results still apply. 
A lower heat flux of 1561 W/m^ was also analyzed from the 
same data run to show the difference between the two heat 
fluxes. All uncertainties are presented as a percentage of 
the calculated parameter, i.e. a relative uncertainty. The 
uncertainty analysis method suggested by Kline and McClintock 

[Ref. 29] was used. For example: If R = R(Xi,X 2 ,X 3 , xj , 

then: 




where : 

6R = uncertainty of the desired dependent variable 
x„ = measured variable 

6x„ = uncertainty in the measured variable 



‘*This analysis is essentially identical to the analysis 
done by Perry [Ref. 2] and therefore his procedure is largely 
reproduced. 
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A. HEAT TRANSFER RATE UNCERTAINTY 



Ir = 1.66 amps 6lr = 0.025 amps 

Vr = 7.76 volts 6Vr = 0.05 volts 

so we have: 



bOh 

Oh 



r/ 6 Vr \2 ^ 




[( VrJ ' 


[ lr)\ 



hQh 

Qh 



r/0.05\2 ^ 


/ 0.025\2] 


\7.76j 


1 1.66 J J 



1.64% 



B. SURFACE AREA UNCERTAINTY 

The following dimensions were taken from 
manufacturer's data sheet, and from actual measurements: 



Do 


= 15.88 


(mm) 


6Do = 


0.1 


(mm) 


L 


= 203.2 


(mm) 


6L = 


0.1 


(mm) 



Ab = Do X L 



6Ab 


[/ 6£)o\2 ^ 




Ab 


H Do j * 


1 ) J 



the 
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6Ab _ 


\i 


' M' + 


/ 0.1 \2] 


Ah 


[(l5, 


, 88/ 


i 203.2/ 



C. WALL SUPERHEAT UNCERTAINTY 

The following temperatures were measured: 
T1 = 19.44 °C 
T2 = 20.03 °C 
T3 = 19.75 °C 
T4 = 19.63 °C 
T5 = 19.46 °C 
T6 = 19.38 °C 
Tavg =19.62 °C 
Standard Deviation = 0.245 °C 
Tsat = 2.20 °C and 6Tsat = 0.01 °C 



Two = Tavg - Qh 



In(^) 

D1 



2 71 'L'kc 



Note: The second term on the right hand side of the above 

equation is the Fourier heat transfer conduction term, and can 
be neglected for the uncertainty analysis since its value is 
less than the standard deviation, so: 

Two = Tavg = 19.62 °C 

AT = Two - Tsat = 19.62 °C - 2.20 °C = 17.42 °C 
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therefore 



6Ar 

AT 



7 6ATVo\2 ^ 


/ bTsat\^' 


(at) ' 


1 AT j J 



6AT 


r/ 0.245\2 ^ 


/ 0. 


o 

to 


Ar 


il7.42) 


(l7 


.42/ J 



1.41% 



D. HEAT FLUX UNCERTAINTY 

//_ {Qh - 2 X Qf) 

^ Ab 

where: Qh = V x I = (194 V) (3.18 A) = 616.92 W 

and 6Qh = (616.92) (0.0164) = 10.12 W 

q" = 59130 W/m^ (measured) 

Ab = 0.010137 

Qh - 2.Qf = (q) (Ab) = (59130) (0.010137) = 599.40 W 
solving for Qf: 

Qf = (616.92 - 599.40)72 = 8.76 W 

assuming the same proportion uncertainty for Qf, we have: 
6Qf = (8.76) (0.0164) = 0.144 W 
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and for the uncertainty: 




hQh \2 ^ / 2hQf \2^/6Ab\2 
Qh - 2Qf] \Qh - 2Qfj \ Ab ) 



1 

2 



5g" 



<7 



// 



[/lO. 




12x0. 


144\: 


\ 599 


•4j 


\ 599 


•4 ) 



+ (0 . 0063)2 2 



1 . 80% 



E. HEAT TRANSFER COEFFICIENT UNCERTAINTY 

h = q"/AT 



therefore, 



h 




I 6AT\2 2 

\ Arj 



=[(0.018)2 +(0.0141)2] 2 



2.29% 



The above series of calculations was also conducted with a 
lower heat flux for the same data run and with a high and low 
heat flux from a TURBO-B data run (D1022FTBC4) with the 
results tabulated in Table 5. 
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TABLE 5. RESULTS OF UNCERTAINTY ANALYSIS 



— 

PARAMETER 


SMOOTH 
1561 W/m" 


SMOOTH 
59130 W/m" 


TURBO- B 
1703 W/m^ 


— 

TURBO-B 
71270 W/m^ 


6Qh/Qh (%) 


8.49 


1.64 


10.47 


1.60 


6Ab/Ab (%) 


0.63 


0.63 


0.63 


0.63 


6AT/AT (%) 


3.55 


1.41 


26.62 


12.44 


6q”/q" (%) 


13.05 


1.80 


9.26 


1.55 


6h/h (%) 


13.52 


2.29 


28.18 


12.54 
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